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CHAPTER  I 


INTRODUCTION 


• 

In  recent  years  delivering  systems  that  perform 
adequately  for  a  specified  period  of  time  in  a 
given  environment  has  become  an  important  goal  for 
both  industry  and  government.  In  the  space 
program,  higher  system  reliability  meant  the 
difference  between  life  and  death.  In  general,  the 
cost  of  maintaining  and/or  repairing  electronic 
equipment  during  the  first  year  of  operation  often 
exceeds  the  purchase  cost,  giving  impetus  to  the 
study  and  development  of  reliability  techniques 
C4:6091 . 

Sound  maintenance  policies  have  their  roots  in 
reliability  theory.  From  reliability  theory,  probability 
models  have  been  developed  for  selecting  optimum  maintenance 
actions.  In  this  context,  maintenance  actions  include 
replacement,  repair,  and  inspection.  The  purpose  of  this 
research  project  is  to  apply  a  theoretical  preventive 
maintenance  model  using  graphical  solution  techniques  to  an 
empirical  data  set  to  determine  optimum  replacement 
intervals  for  synthetic  and  petroleum  motor  oils. 


Preventive  maintenance  is  defined  as  follows: 

Equipment  maintenance  actions  performed  on  a 
periodic  basis,  according  to  a  specific  set  of 
instructions  and  a  predetermined  time  schedule. 
The  objective  is  to  protect  equipment  capability 
and  investment  by  removing  the  causes  of  failure 
and  making  adjustments  to  compensate  for  normal 
wear  before  failure  occcurs  E 12:5351. 
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In  other  words,  preventive  maintenance  is  considered  as 
those  actions  designed  to  preclude  -failure  and  enhance 
system,  subsystem,  or  component  performance.  On  the  other 
hand,  corrective  maintenance  can  be  considered  as  those 
actions  designed  to  restore  a  system,  subsystem,  or 
component  and  to  enhance  performance  after  failure  has 
occurred.  Further,  the  authors  consider  an  optimum 
maintenance  strategy  as  one  that  maximizes  or  minimizes  the 
objective  function  of  the  maintenance  policy,  such  as 
availability  or  cost.  Thus,  an  optimum  maintenance  strategy 
could  be  either  preventive  or  corrective. 

Preventive  maintenance  policies  can  be  categorized 
according  to  uncertainty  about  failure.  This  project 
addresses  three  types  of  uncertainty.  The  first  two  types 
are  defined  here,  and  the  third  type  will  be  discussed  as 
part  of  the  problem  statement. 

1.  Type  one  uncertainty  exists  when  there  is  dt-ubt 
regarding  the  component's  life  (i.e.,  when  the  next 
failure  will  occur).  This  type  of  uncertainty  is 
present  for  all  components  demonstrating  stochastic 
failure  patterns. 

2.  Type  two  uncertainty  exists  when  there  is  doubt 
regarding  the  present  condition  of  the  component 
<i.e.,  whether  the  component  is  good  or  failed). 

The  presence  of  this  type  uncertainty  depends  on  the 
nature  of  the  component. 

In  cases  where  only  type  one  uncertainty  is  present,  and  the 
component's  condition,  good  or  failed,  is  Known  at  any  point 
in  time,  either  a  planned  replacement  or  repair  policy  may 
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be  optimal.  In  casts  where  both  type  one  and  type  two 
uncertainty  are  present,  and  the  component's  condition  is 
not  known  with  certainty,  either  a  planned  inspection  or 
what  Talbott  (11)  has  termed  the  ‘blind  replacement"  policy 
may  be  the  optimum  preventive  maintenance  strategy. 

However,  depending  upon  the  component's  -failure  pattern  and 
the  pay-o-ffs  associated  with  maintenance  actions  versus 
failure,  prevent i ve  maintenance  may  make  no  sense.  The 
optimum  policy  may  be  to  take  no  action  until  failure  and 
then  perform  correct i ve  maintenance. 


Problem  Statement 

While  theoretical  probability  models  exist  for 
selecting  optimum  maintenance  policies,  these  techniques  are 
often  not  adaptable  to  real-world  preventive  maintenance 
policies;  thus,  they  receive  little  recognition  in  Air  Force 
maintenance  policies.  A  significant  shortcoming  of  most 
probability  models  is  that  they  require  rather  explicit 
knowledge  of  the  component's  underlying  failure 
di str i but i on . 

All  statistical  procedures  require  an  unbiased 
sample  of  a  population.  To  the  extent  that  the  sample  does 
not  exactly  model  the  population,  a  potential  for  error 
exists  when  drawing  conclusions  based  on  sample  data.  When 
this  sample  is  component  failure  data,  the  error  is  a  direct 
result  of  the  type  one  uncertainty  present  in  all  stochastic 


processes.  Additionally,  if  in  order  to  apply  a  particular 
methodology,  one  must  be  able  to  describe  the  sample  data  in 
terms  of  some  theoretical  distribution,  then  further 
potential  for  error  exists  when  drawing  conclusions  based  on 
the  theoretical  distribution.  This  second  source  of  error 
is  a  direct  result  of  the  extent  to  wh.ch  the  theoretical 
distribution  does  not  exactly  describe  the  sample  data. 

Even  when  ample  failure  data  is  available,  if  failure 
patterns  do  not  fit  theoretical  failure  distributions,  most 
probability  models  cannot  reasonably  be  applied.  Hence,  the 
problem  becomes  one  of  development  and  application  of  sound 
maintenance  strategies  that  addresses  a  third  type  of 
uncertainty:  uncertainty  about  the  underlying  failure 

distribution  of  the  component.  Thus,  type  three  uncertainty 
exists  when  there  is  doubt  regarding  the  failure 
distribution  of  the  component. 

S&SBJL 

The  authors  believe  an  appropriate  solution  to  this 
problem  is  to  apply  graphical  solution  techniques  to 
preventive  maintenance  models.  If  a  particular  methodology 
does  not  require  knowledge  of  the  theoretical  distribution, 
error  resulting  from  type  three  uncertainty  is  eliminated. 
Graphical  techniques  do  not  require  fitting  a  theoretical 
distribution  to  the  component's  failure  distribution;  hence, 
they  eliminate  error  associated  with  type  three  uncertainty. 
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The  particular  graphical  technique  we  will  use  is  the  total 
time  on  test  <TTT)  transform  which  Bergman  (3:467-469)  has 
already  adapted  for  a  preventive  maintenance  policy  Known  as 
the  age  replacement  model. 

According  to  Talbott  <11>,  techniques  for  adapting 
the  TTT  transform  to  other  preventive  maintenance  policies 
have  yet  to  appear  in  the  literature  but  are  currently  under 
investigation.  Additionally,  application  of  these  solutions 
to  real-world  situations  using  an  empirical  data  set  remains 
to  be  demonstrated.  One  particular  maintenance  model, 
termed  the  blind  replacement  model,  proposes  replacement  of 
components  at  planned  intervals  without  Knowledge  of 
component  condition  (4:190-192).  Talbott  (10:16)  is 
developing  a  graphical  solution  technique  for  this 
maintenance  model.  Uie  believe  application  of  this  technique 
to  an  empirical  data  set  will  demonstrate  the  validity  of 
using  graphical  solutions  for  developing  real-world 
maintenance  strategies. 

A  disadvantage  of  TTT  analysis  is  that  it  probably 
will  not  identify  the  truly  optimal  solution.  The  nature  of 
this  technique  is  to  select  from  the  observed  failure  data, 
the  data  point  that  is  most  optimal  among  all  the  data 
points  for  the  preventive  maintenance  tasK .  Thus,  unless 
there  is  an  observed  failure  at  the  optimal  preventive 
maintenance  interval,  the  model  will  not  select  the  optimal 
interval,  but  rather,  the  data  point  closest  to  it.  In 
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other  words,  the  model  considers  the  empirics!  date  set  as  a 
finite  set  of  potential  replacement  intervals  from  which  the 
optimal  interval  can  be  selected.  We  will  elaborate  on  this 
disadvantage  as  we  discuss  the  TTT  methodology.  However,  in 
the  interest  of  simplification,  we  wilt  continue  to  address 
the  interval  selected  by  TTT  analysis  as  optimal  even  though 
it  may  not  be  truly  optimal. 

Our  empirical  data  base  for  this  investigation  is 
synthetic  and  petroleum  lubricant-lifetime  data  made 
available  by  the  Technical  Support  Center,  Pensacola  Naval 
Air  Station,  FL.  This  data  base  is  part  of  a  joint-service 
test  of  DOD  vehicles  to  determine  the  merits  of  both  types 
of  lubricants.  Application  of  the  blind  replacement  model 
using  graphical  solutions  should  yield  optimal  oil 
replacement  intervals  for  synthetic  versus  petroleum  motor 
oils. 

Research  Question 

Can  the  blind  replacement  mode)  be  applied  to  an 


empirical  data  set  using  TTT  transform  graphical  solution 
techniques  to  determine  optimal  replacement  intervals  for 
synthetic  and  petroleum  motor  oils? 


CHAPTER  II 


LITERATURE  REVIEW 

This  chapter  presents  a  brief  look  at  the 
development  of  optimal  maintenance  policies  based  on 

reliability  theory.  Barlow  and  Proschan  <2i84>  categorize 
maintenance  policies  as  replacement  policies  and  inspection 
policies.  The  topics  covered  in  this  chapter  are 
reliability  theory,  maintenance  policies  in  general, 
replacement  models,  inspection  models,  and  the  blind 
replacement  model.  Discussion  of  the  TTT  transform 
graphical  technique  and  its  application  to  maintenance 
policies  is  reserved  for  Chapter  Three,  Methodology. 

BtlUfailliz  Ih*.acr 

Barlow  and  Proschan  <2*6>  define  reliability  as  "the 
probability  of  a  device  performing  its  purpose  adequately 
for  the  period  of  time  intended  under  the  operating 
conditions  encountered."  Similarly,  Hillier  and  Lieberman 
<4*394)  state  that 

reliability,  R<t>,  is  the  probability  that  a  device 
performs  adequately  over  the  interval  <0,t>.  In 
general,  it  is  assumed  that  unless  repair  or 
replacement  occurs  adequate  performance  at  time  t 
implies  adequate  performance  during  the  interval 
<0,t>.  The  device  under  consideration  may  be  an 
entire  system,  a  subsystem,  or  a  component. 

Some  probability  models  of  reliability  consider 
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devices  to  bf  in  on*  of  two  conditions*  good  or  failed.  For 
an  individual  component,  this  condition  can  be  denoted  by  a 
random  variable  Xj f  such  that  Xj  *1  i f  the  i  component 
is  good  and  Xj  *  0  if  the  i  component  is  failed.  Thus,  Xj 
is  a  binary  random  variable  defined  as; 

!1,  if  device  i  does  not  fail  during  <0,t) 

0,  if  device  i  does  fail  during  (0,t> 

Therefore,  the  probability  of  a  device  not  failing  during 
the  interval  <0,t)  can  be  expressed  as  P<X.(t>s*l>  or  P<X»1> 
<4*595) . 

The  common  probability  functions  used  in  reliability 
theory  treat  system  condition  as  a  function  of  time.  The 
function  f<t>  represents  the  failure  probability  density 
function  <pdf) — that  is,  the  probability  of  failure  at  time 
t;  and  F<t)  represents  the  cumulative  failure  pdf — that  is, 
the  probability  of  failure  within  the  interval  <0,t).  The 
function  R<t)  represents  the  survival  distribution  or 
reliability  distribution  and  is  the  probability  of  survival 
throughout  the  interval  <0,t>;  hence, 

R<  t)  -  P<X«1>  -  1 -F<  t >  (2.1) 

The  hazard  rate,  r<t>,  sometimes  referred  to  as  the 
instantaneous  failure  rate,  is  the  conditional  probability 
that  a  device  surviving  to  time  t  wi 1 1  fail  at  time  t. 
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Tharafora,  it  is  tha  ratio  o-f  tha  probability  of  failure  at 
tima  t  to  tha  probability  of  survival  throughout  tha 
intarval  <0,t)f  hanca, 

r<  t>  -  f <  t>/R<  t)  <2. 2) 

For  many  davicas,  tha  hazard  rata  tands  to  incraasa  as  tha 
davica  gats  oldar  dua  to  da tar i or at i on  assoc iatad  with  aga. 
Oavicas  having  hazard  ratas  that  ramain  constant  or  incraasa 
with  aga  ara  said  to  hava  an  incraasing  failura  rata  (IFR). 
For  soma  davicas,  tha  hazard  rata  tands  to  dacraasa  as  tha 
davica  gats  oldar  dua  to  a  possibla  'work  hardaning" 
phanomanon  that  occurs  during  aarly  lifa.  Davicas  having  a 
hazard  rata  that  ramains  constant  or  dacraasas  with  aga  ara 
said  to  hava  a  dacraasing  failura  rata  <DFR).  Davicas 
having  an  axponantial  failura  distribution  hava  a  random  or 
poisson  failura  rata  that  is  constant  throughout  tima; 
hanca,  tha  axponantial  distribution  is  both  IFR  and  DFR  and 
raprasants  a  natural  boundary  on  tha  survival  probabi 1 i t i as 
of  tha  IFR  and  DFR  distributions  (4i605-606>. 

Sinca  R(t>  is  a  continuous  function  which  cannot 
assuma  nagativa  valuas,  tha  maan  valua  of  X  (axpactad 
condition  of  tha  davica)  is  raprasantad  by  tha  following 
axprassion  (5:57) i 

oo 

E<x)  -  (  R(x)dx  (2.3) 

o' 
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Furthermore,  -failure  distributions  represent 
probability  models  of  the  length  of  life  of  e  component  or 
system.  Typical  failure  density  functions  popular  in  the 
1 i terature  are  the  exponential ,  the  gamma,  the  Wei  bull,  the 
normal,  and  the  lognormal.  The  exponential  distribution  has 
limited  utility  <a!beit  wide  application)  because  of  its 
underlying  property  that  previous  use  does  not  impact  future 
life.  In  other  words,  a  new  device  has  the  same  probability 
of  failing  within  the  next  time  interval  as  an  old  device. 
For  devices  having  exponential  failure  distributions,  the 
future  life  distribution  of  the  devi-  *  remains  unchanged  as 
long  as  the  device  has  not  failed  yet.  Therefore,  devices 
having  exponential  failure  distributions  essentially  do  not 
age  <2i9-lS). 

Maintenance  Policies 

For  some  devices,  failure  during  operation  is 
undesirable  because  it  is  too  costly  or  dangerous.  In  these 
situations,  it  may  be  wise  to  employ  a  maintenance  policy 
that  incorporates  planned  repl acemen t/i nspec t i on  of  good 
devices.  Barlow  and  Proschan  <2i84>  define  an  optimum 
maintenance  policy  as  one  'that  minimizes  total  cost, 
maximizes  availability,  or  in  general  attains  the  best  value 
of  the  prescribed  objective  function." 


Central  to  the  development  of  maintenance  policies 
is  renewal  theory.  A  renewal  process  may  be  thought  of  as 


any  process  that  returns  a  device  to  a  good-as-new 
condition |  this  includes  both  replacement  and  repair  to  a 
good-as-new  condition.  The  renewal  reward  theorem  states 
that  the  average  long-term  reward  is  equal  to  the  expected 
reward  in  a  renewal  cycle  divided  by  the  expected  cycle 
length.  Mathematically,  this  is  expressed  as: 

E<x/cycl e) 

E<x)  - -  (2.4) 

E(cyc1 e  length) 

Here  X  represents  some  -form  o-f  reward  or  payoff  which  can  be 
positive,  such  as  availability,  or  negative,  such  as 
dawn-time  or  cost.  The  objective  function  then  becomes 
maximization  of  a  positive  long-term  reward  or  minimization 
of  a  negative  long-term  reward  <7iS3>. 

In  this  research  project,  we  will  use  the  term 
cost  to  signify  any  negative  reward,  and  prof i t  for  any 
positive  reward.  Also,  the  terms  renewal  and  replacement 
are  used  interchangeably,  such  that  replacement  refers  to 
either  physical  exchange  of  a  device  or  repair  to  a  good- 
as-new  cond i t i on . 

Bujimmiil  Motif  If 

Preventive  maintenance  policies  are  categorized  by 
Barlow  and  Proschan  <2s84)  as  either  replacement  type  or 
inspection  type  policies.  Replacement  policies  attempt  to 
balance  the  cost  of  failure  occurring  during  operation  of 
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the  device  against  th«  cost  of  planned  replacement.  The 
following  assumptions  make  application  of  a  replacement 
policy  reasonable!  1)  cost  of  replacing  a  good  device  is 
less  than  cost  of  the  device/s  failing!  2)  detection  and 
replacement  of  failed  devices  is  instantaneous!  and  3)  the 
probabi 1 i ty  of  failure  is  increasing  with  age.  In  the 
Barlow  and  Proschan  models,  C<t)  represents  the  expected 
long-term  cost  incurred  during  the  interval  (0,t>!  C, 

represents  the  cost  of  device  failure  to  include  cost  of 
replacing  the  failed  device!  and  C2  denotes  the  cost  of 
replacing  a  good  device.  Nj(t)  represents  the  expected 
number  of  failures  occurring  during  the  interval  (0,t>, 
while  N2<*>  represents  the  expected  number  of  replacements 
of  good  devices.  The  basic  replacement  model  can  then  be 
expressed  asi 


C<t>  -  CjN^t)  ♦ 


(2.5) 


Here,  the  objective  is  to  minimize  C(t)  for  a  finite  time 
span  or  minimize  |jj^C<t)/t  for  an  infinite  time  span 


(2184-85) . 


The  most  common  replacement  model  is  the  age 
replacement  model  which  proposes  replacing  devices  at 
failure  or  at  age  T,  whichever  occurs  first  (where  T  is 
assumed  to  be  a  constant).  Barlow  and  Proschan  show  that  if 
the  underlying  failure  distribution  is  continuous,  then 
determining  an  optimum  age  of  replacement  requires  exact 


12 


Knowledge  of  the  failure  distribution,  linfor tunatel y ,  in 
actual  practice,  the  -failure  distribution,  F,  may  not  be 
Known.  However,  if  the  failure  distribution,  F,  is  Known 
exactly,  the  optimum  age  of  replacement  for  an  infinite  time 
span  is  that  value  of  T  which  yields  the  minimum  value  for 
C(T>  in  the  following  equation  (2i85-86). 


C.F(T>  ♦  C,R(T) 

C<T)  «  — - - - - 

T  (2.6) 

(  R(x>dx 

0 


This  equation  can  be  derived  mathematically  from 
equation  2.S  and  is  equivalent  to  the  mathematical 
expression  of  the  renewal  reward  theorem  (equation  2.4) 
where i 

C(T)  *  expected  long-term  cost  at  time  T 

C,F(T>  ■  failure  cost  multiplied  by  the  probabi 1 i ty 

1  of  failure  within  the  interval  <0,T> 

C_R(T)  *  planned  replacement  cost  multiplied  by  the 

2  probability  of  survival  throughout  the 
interval  <0,T) 

T 

(  R<x)dx  *  the  expected  value  of  the  cumulative  survi- 
0J  val  pdf  where  R(x)  is  bounded  by  0  and  T 

(equation  2.3) 

Hence,  C,F(T)  ♦  C9R(T)  ■  E(cost  per  cycle),  and 
T  1  * 

^  R(x)dx  ■  E(cycle  length)  (2i8S-86). 

0 

Regarding  finite  time  spans,  Barlow  and  Proschan 
(2»?3>  show  that  if  F  is  Known  and  continuous,  then  there 
exists  an  age  of  replacement  for  any  finite  time  span  (0,T){ 
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howvvtr ,  no  general  formula  exists  for  determining  an 
optimum  age  replacement  policy. 


A  second  type  of  replacement  policy  is  block 
replacement,  which  proposes  replacing  devices  at  calendar 
times  kT  <k  -  1,2...),  independent  of  component  failure 
history)  devices  are  also  replaced  at  failure.  Block 
replacement  policies  have  an  advantage  over  age  replacement 
in  that  they  eliminate  the  administrative  burden  of  tracking 
device  age.  However,  the  disadvantage  of  a  block 
replacement  policy  is  that  it  throws  away  useful  life.  In 
comparing  block  and  age  replacement  policies,  Barlow  and 
Proschan  show  that  with  block  replacement,  more  unfailed 
devices  are  removed,  more  useful  life  is  wasted,  and  the 
total  number  of  both  failed  and  good  device  removals  is 
greater  than  with  age  replacement)  however,  the  expected 
number  of  failures  is  less  for  block  replacement  policies 
<2* 67,93) . 

iDlBMtlOP  flarftl* 

As  stated  earlier,  a  basic  assumption  of  replacement 
policies  is  that  detection  and  replacement  of  a  failed 
device  are  instantaneous.  Thus,  replacement  models  assume 
that  the  condition  of  the  device  (failed  or  good)  is  known 
at  any  point  in  time.  When  type  two  uncertainty  exists 
<i.e.,  the  condition  of  the  device  cannot  be  determined 
instantaneously),  replacement  policies  cannot  logically  be 
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applied;  and  inspection  policies,  which  make  up  the  second 
major  category  of  maintenance  policies,  are  more  suited  for 
these  i terns  ( 1 1 > . 

Inspection  policies  attempt  to  balance  the  cost  of 
failure  occurring  during  operation  of  the  device  against  the 
cost  of  planned  inspections.  The  following  assumptions  make 
application  of  an  inspection  policy  reasonable:  1>  the 
condition  of  the  device  can  be  determined  only  through 
inspection;  2)  inspection  is  non-destructive;  3)  the  device 
cannot  fail  during  inspection;  4)  inspection  involves 
a  fixed  cost,  j  *nd  the  time  elapsed  between  system 
failure  and  discovery  of  failure  at  the  next  inspection 
costs  C2  per  unit  of  time.  As  with  the  age  replacement 
model,  Barlow  and  Proschan  apply  a  mathematical  model  to 
describe  the  relationship  between  the  expected  long-term 
cost  C(t)  and  Cj  and  C2*  They  define  an  optimum  inspection 

procedure  as  a  'specification  of  specific  checks  Xj  (. 

X_  <  X<a**»*  'for  which  C(T)  is  minimized;  and  they  show 
2  d 

that  if  the  failure  distribution,  F,  is  known  and  continuous 
with  a  finite  mean,  then  an  optimum  inspection  schedule  must 
exist  (2:107-108). 

Blind  Replacement  Model 

Radner  and  Jorgenson  (6:184)  categorize  maintenance 
models  according  to  sources  of  uncertainty.  One  source  of 
uncertainty  occurring  with  all  maintenance  policies 
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developed  -for  stochastically  -failing  devices  is  predicting 
when  a  -failure  will  occur  (type  one  uncertainty).  A  second 
source  of  uncertainty  is  knowing  the  condition  of  the  device 
(good  or  failed)  at  any  point  in  time  without  inspection 
(type  two  uncertainty).  Barlow  and  Proschan  apply 
replacement  models  to  cases  where  only  type  one  uncertainty 
is  present  and  apply  inspection  models  to  cases  where  both 
type  one  and  type  two  uncertainty  are  present.  Similarly, 
Radner  and  Jorgenson  propose  a  category  of  preparedness 
models  for  situations  having  both  types  of  uncertainty. 
Preparedness  models  include  but  are  not  restricted  to  the 
typical  inspection  type  models  addressed  by  Barlow  and 
Proschan.  For  example,  one  assumption  Barlow  and  Proschan 
make  regarding  inspection  models  is  that  inspection  is 
nondestructive.  However,  in  some  situations,  inspection  may 
cause  the  device  to  fail.  This  is  true  of  any  device  that 
must  be  expended  in  order  to  determine  if  the  device  is 
still  good,  such  as  munitions  which  must  be  fired  to 
determine  if  they  still  work  but  clearly  will  not  work  again 
after  they  have  been  fired.  For  devices  falling  in  this 
category,  planned  inspection  makes  no  sense,  and  the  optimum 
policy  would  be  to  replace  at  a  planned  interval.  However, 
these  devices  do  not  meet  the  requisites  of  Barlow  and 
Proschan's  reel acement  policies  since  both  type  one  and 
type  two  uncertainty  are  present. 

Radner  and  Jorgenson  developed  a  preparedness  model 
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-for  devices  hewing  an  arbitrary  distribution  of  times  to 
failure,  which  Talbott  (11)  has  termed  the  'blind 
replacement*  model.  Blind  replacement  proposes  renewing  a 
device  at  age  or  time  T  without  knowledge  of  the  condition. 
Assumptions  for  the  blind  replacement  model  are  !)•  the  cost 
of  replacement  is  less  than  cost  of  inspection  or  inspection 
is  an  unreasonable  option,  and  2)  the  probability  of  failure 
is  increasing  with  age.  If  A(T>  represents  the  average 
long-term  profit  (availability,  in  this  case);  T,  the  time 
between  replacements;  and  k,  the  time  to  replace,  then  A(T) 
can  be  expressed  in  terms  of  reward  renewal  theorem 
(equation  2.4)  as: 


(  R(x)dx 

0  (2.7) 

A(T)  -  - 

T  +  K 

Here,  the  numerator  of  the  above  equation  represents  the 
mean  life  in  a  renewal  cycle;  and  the  denomonator  represents 
the  time  between  replacements  plus  time  to  replace,  which  is 
cycle  length.  The  objective  function  then  is  to  maximize 
A(T>  (6:190-192). 

Summary 

Two  important  requisites  exist  for  all  the 
maintenance  policies  presented  thus  far.  First,  the  failure 
distribution  must  be  known;  and  second,  it  must  be  IFR,  but 
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not  exponent i *1 .  Although  Barlow  and  Proschan  provide 
several  mathematical  proofs  of  the  requisite  for  an  IFR 
failure  distribution,  it  is  intuitive  that  devices  having 
purely  exponential  or  DFR  distributions  are  not  reasonable 
candidates  for  preventive  maintenance.  Since  devices  having 
purely  exponential  failure  distributions  essentially  do  not 
age  and  have  a  constant  hazard  function  throughout  time, 
then  clearly,  replacing  a  good  device  will  not  impact  the 
probability  of  failure  within  the  next  time  interval;  in 
other  words,  the  new  device  would  have  the  same  probability 
of  failing  as  the  old.  For  devices  having  non-exponential, 
OFR  failure  distributions,  planned  replacement  of  a  good 
device  actually  increases  the  probability  of  failure  within 
the  next  time  interval  since  the  new  device  would  have  a 
greater  probability  of  failure  than  the  old.  Therefore,  an 
underlying  assumption  of  all  preventive  maintenance  policies 
discussed  in  this  project  is  that  their  application  is 
reasonable  only  if  the  device  has  a  non-exponential,  IFR 
failure  distribution. 

A  major  drawback  of  all  the  maintenance  models 
presented  thus  far  is  that  determining  an  optimum 
rep 1 acement/ inspect i on  strategy  requires  an  exact 
description  of  the  underlying  failure  distribution,  F,  in 
terms  of  some  theoretical  failure  distribution  such  as  a 
gamma  or  Uieibull.  As  stated  earlier,  in  actual  practice, 
the  failure  distribution  may  not  be  known  or  may  not 
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reasonably  fit  a  theoretical  distribution.  This  fact  alone 
makes  application  of  theoretical  maintenance  policies 
difficult  for  most  real-world  preventive  maintenance 
situations.  Thus,  these  models  could  be  made  more  useful  if 
they  could  be  adapted  to  data  analysis  techniques  employing 
empirical  distributions.  One  such  technique  is  a  graphical 
technique  using  total  time  on  test  transform,  which  is  the 
subject  of  the  next  chapter. 


CHAPTER  III 


METHODOLOGY 

This  chapter  presents  the  total  time  on  test 
methodology  for  analyzing  observational  data  to  determine 
optimal  maintenance  intervals.  This  technique  was  described 
by  Barlow  and  Campo  <1:451)  as  an  effective  means  of 
analyzing  failure  data.  In  1977,  Bergman  (3:467-469) 
adapted  Barlow  and  Proschan's  age  replacement  model  to  TTT 
analysis  techniques.  Also,  Talbott  (10:16)  is  currently 
adapting  Radner  and  Jorgenson's  blind  replacement  model  for 
TTT  analysis.  In  describing  the  development  of  the  TTT 
transform  technique,  this  chapter  will  cover  the  TTT  process 
in  general,  the  Bergman  model,  the  Talbott  model,  and  data 
appl i cat  ion. 

Total _Time  on  Test  Process 

The  total  time  on  test  (TTT)  process  was  presented 
by  Barlow  and  Campo  in  1975  as  a  method  of  analyzing  data 
that  is  useful  for  describing  failure  probability.  An 
explanation  of  this  process  follows: 

*(j>  is  a  value  from  an  ordered  set  of  life 
observations  from  distribution  F  such  that  t^jj  £  t(2>  1 

*<3)  i  •**  i  *<n)‘  Since  t^j^  represents  the  life  of  the 
first  device  that  failed,  and  t<2)  th#  life  of  the  second, 
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etc.,  then  it  -follows  that  n  devices  si)  lived  t^^  units 
o-f  time,  n-1  devices  lived  t(2>  units  o f  time,  snd  only  one 
device  lived  t<n)  units  o-f  time.  Thus  t(j)  represents  the 
li-fe  o-f  the  i  device  in  the  ordered  set  o-f  n  lifetimes, 
where  t(0)*0. 

2.  T(  is  the  totsl  time  on  test  statistic  snd 
represents  the  life  generated  by  the  n  devices  in  the 
interval  (t(Q),t<|>>»  where  a  portion  of  the  n  devices 
survived  throughout  the  interval  and  the  remainder 
failed  during  the  interval.  Where  T^“0 ,  then  Tj  is 
defined  mathematically  as, 

i 

ti  "  YL <n  -  ■*  * >>  (tu>  -  ‘(j-D1  <3-n 

j-i 

3.  Uj  is  the  scaled  total  time  on  test  statistic 
and  is  equal  to  the  ratio  of  T}  to  •  Thus,  Uj  represents 
the  proportion  of  the  total  life  generated  by  the  n  devices 
during  the  interval  (t<0)  ’*<  i  >>  ' 

4.  i/n  Is  the  ratio  of  i  to  n  and  represents  the 
proportion  of  the  devices  that  have  failed  by  point  in  time 
t(i).  Thus,  it  is  approximately  equal  to  the  probability 
that  an  item  fails  within  the  interval  *t<o>**<i>*' 

9.  The  TTT  plot  is  the  graphical  representation  of 
the  TTT  process  where  U(  is  plotted  against  i/n.  Therefore, 
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each  point  on  the  TTT  plot  relates  the  proportion  of  the 
total  Ufa  ganaratad  by  tima  t<n  to  tha  probability  that  a 
davica  fails  by  that  tima.  Necessarily  than,  at  tha  final 
point  where  i«n,  Uj  ■  i/n  ■  1  <li452-437}  11). 

Barlow  and  Campo  <li453>  furthar  show  that  if  tha 
data  points  on  tha  TTT  plot  ara  joinad  by  lina  sagmants, 
whan  n  is  sufficiantly  larga,  this  function  will  approxfmata 
a  curva.  Tha  curve  will  ba  concave  for  1FR  distributions,  a 
straight  lina  for  puraly  axponantial,  and  convax  for  DFR 
distributions.  Our  axamplas  of  thasa  functions  ara  givan  in 
appandicas  A  through  C  of  this  report. 

Tha  Bergman  Modal 

Bergman  made  a  major  breakthrough  in  tha  application 
of  probability  based  maintenance  strategies  to  real-world 
scenarios  by  adapting  Barlow  and  Proschan's  age  replacement 
modal  to  TTT  transform  analysis.  This  method  provides  a 
useful  means  of  estimating  tha  optimum  age  replacement 
interval  whan  only  observational  data  is  available,  and  tha 
device's  underlying  failure  distribution  is  not  known  with 
certainty  or  cannot  ba  fitted  to  an  empirical  distribution. 
Tha  modal  assumes  a  cost  C  is  associated  with  component 
replacement  and  an  additional  cost  C+K  is  associated  with 
failure.  If  we  let  K*1 ,  than  from  aquation  2.6,  tha  optimum 
age  of  replacement  is  that  value  of  T  which  yields  tha 
minimum  value  for  C(T)  in  tha  following  aquations 
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(CMXF(T>>  ♦  C(R<T)> 

C<T)  -  - 

T  (3.2) 

(  R< x >dx 


Sine*  R<T>  ■»  1-F(T>,  the  numerator  of  the  above  equation  can 
be  reduced  mathematically  to  C+F(T> .  Also  since  Fn<T)  is  an 
empirical  distribution  o-f  n  observations  that  approaches  the 
theoretical  F<T)  as  n  approaches  infinity,  then 
FCT)  ■  1 im  F  (T>  .  Further,  since  represents  one 

observation  -from  that  empirical  distribution,  then  (using 


to  represent  "is  approximately  equal  to") 

F(T>  ■*  F  (T>  F  (t>..)  when  n  is  su-f-f  ici  ently  large.  And 
n  ~  n  '  i  * 

necessarily,  C<T)  &  cn<T)  “  Cn<t<i>>*  Hence* 


C(T> 


C  ♦  F<T) 
T 

(  R(x)dx 

o' 


<*&'ix>dx 
0 '  n 


Cn(t<  i  >  * 


(3.3) 


Since  T|  represents  the  cumulative  li-fe  generated  by  the 

n  components  over  the  interval  **(□>**< i>*  *nd  since 
,  *<  i  > 

S  R  (x)dx  represents  the  expected  value  of  the  cumulative 

O'  n 

survival  pdf  over  that  same  interval,  then  for  the  empiri- 

t  /  j  \ 

cal  data  set,  $  RR(x)dx  -T,/n.  Also,  since  t(j)  repre¬ 
sents  one  data  point  in  the  empirical  cumulative  failure  pdf, 
Fn(T>,  then  necessarily,  the  probability  of 

failure  within  the  interval  *nd  •*  therefore 

approx imatel y  equal  to  i/n.  Thus  from  equation  3.3,  it  can 
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be  shown  that, 


C  ♦  i/n 

-  <3.4> 

< i/n>  <T  >  <T,/T  ) 


Note  that  th*  -first  part  o-f  th*  expression  on  th*  right 
hand  side  is  a  constant  fixed  by  th*  sample,  and  the 
second  part  is  the  only  part  that  is  variable.  Sine*  the 
objective  is  to  minimize  long-term  costs,  this  is  equivalent 
to  maximizing  th*  reciprocal  o-f  Cn<t(.}).  Therefore  the 
objective  becomes  maximizing  the  reciprocal  of  the  variable 
portion  of  th*  right  hand  side  of  equation  3.6.  Hence, 

min  C<T>  — >  max - 1 -  (3.7) 

C  +  i/n 


For  a  particular  value  of  i,  plotting  -C  on  th*  X-axis  of 
th*  TTT  plot  yields  a  distance  from  -C  to  i/n  that  equals  C 
♦  i/n.  If  a  line  is  drawn  from  -C  on  th*  X-axis  to  th*  TTT 
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data  point  correspond! ng  to  that  value  of  i,  that  line  is 
defined  by  its  slope  which  is  equal  to  U(/(C  ♦  i/n).  For 
each  line  that  can  be  drawn  in  this  fashion,  at  least  one 
line  will  be  tangent  to  the  function  and  will  have  a  slope 
that  is  greater  than  or  equal  to  the  slope  of  any  other 
line  that  can  be  drawn.  Thus,  that  line  defines  the  value 
of  i  (call  it  i«>  that  maximizes  Uj/(C  ♦  i/n)  for  a  given 
value  of  C.  The  optimal  replacement  interval  then  is  the 
value  of  t(j)  defined  by  i»  <3i447-449j  8!?-10>. 

As  stated  earlier,  the  disadvantage  of  TTT  analysis 
is  that  it  probably  will  not  identify  the  truly  optimal 
solution.  Since  the  technique  selects  the  value  of  t(j) 
defined  by  the  point  of  tangency,  it  can  consider  only  the 
observed  lifetimes  as  candidates  for  the  replacement 
interval.  If  the  truly  optimal  solution  is  not  a  value  of 
t^jj,  the  point  of  tangency  will  be  the  ttj)  closest  to  the 
truly  optimal  interval. 

Additionally,  if  the  distribution  is  purely 
exponential,  the  TTT  function  will  be  a  straight  line,  and 
if  it  is  DFR,  the  function  will  be  convex.  In  either  of 
these  two  cases,  the  point  of  tangency  will  be  the  last  data 
point  on  the  TTT  plot.  Hence,  the  model  should  be 
interpreted  for  purely  exponential  and  DFR  distributions  as 
having  an  optimal  strategy  of  replacing  only  at  failure. 

Our  examples  of  IFR,  purely  exponential  and  DFR 
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distributions  ara  givan  in  Abtachmants  A  through  C  of  this 
report . 

Tha  Talbott.  Hods  1 

From  Bargman's  insight,  Talbott  <10il6)  is  adapting 
Radnar  and  Jorganson's  blind  raplacamant  modal  to  TTT 
analysis  tachniquas.  Using  a  procass  similar  to  Bargman's, 
Talbott  shows  that  as  n  approach#*  infinity,  tha  ampirical 
distribution  approachas  tha  thaoratical  distribution.  Thus, 
from  aquation  2.7,  ha  shows  that 
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of  both  types  of  lubricants.  Application  of  the  blind 
replacement  modal  using  graphical  solutions  should  yield 
optimal  oil  replacement  intervals  for  synthetic  versus 
petroleum  motor  oils. 
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CHAPTER  IV 


DATA  ANALYSIS 

This  chapter  examines  the  procedures  used  to  collect 
and  analyze  the  oils  comprising  the  data  base  -for  this 
research.  It  also  discusses  the  method  -for  determining  K, 
the  time  required  to  replace  -failed  oil. 

The  original  motivation  -for  the  data  collection 
e-f-fort  was  to  determine  i -f  there  were  any  advantages  to 
using  synthetic  lubricants  in  DoD  vehicles.  Air  Force 
Systems  Command  (AFSC)  was  tasked  to  per-form  a  -field 
evaluation  to  determine  i-f  using  synthetic  engine  oils 
provided  improved  performance,  greater  reliability,  and 
lower  vehicle  operating  costs  <9sl>. 

Data  Col  lection 

Oil  sample  data  -for  the  field  evaluation  was 
collected  over  a  three-year  period  <1979-1981)  through  the 
Management  and  Equ  ment  Evaluation  Program  <MEEP>  Project 
H79-1C.  The  test  included  approx imate 1 y  450  general  purpose 
vehicles  from  twelve  Air  Force  bases  throughout  the  country. 
The  large  number  of  vehicles  and  the  wide  cro .s-sec t i on  of 
bases  should  provide  a  relatively  unbiased  simple  of  Air 
Force  vehicles.  The  lubricants  selected  for  i'.e  program 
were  MIL-L-46152  Qualified  Products  List  <QPL>  products. 
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Three  different  types  of  lubricants  were  evaluated: 


Stauffer's  synthetic  oil,  a  standard  stock  listed  mineral 
based  oil,  and  CONOCO' s  synthetic  oil.  To  insure 
consistency  in  determining  oil  failure,  oil  samples  were 
taken  at  regular  intervals  and  sent  to  the  Joint  Oil 
Analysis  Program  <JOAP)  Laboratory  at  the  Naval  Air  Station, 
Pensacola,  Florida  for  testing.  Oil  filters  were  changed 
each  6,000  miles  or  12  months  of  operation,  whichever  came 
first.  The  personnel  at  the  JOAP  Laboratory  made 
recommendations  based  on  their  findings  whether  or  not  to 
change  oil/filter  at  more  frequent  intervals  or  to  continue 
with  routine  sampling  <9:i-6). 

Oil  was  considered  failed  when  certain  established 
baseline  factors  were  exceeded.  Factors  included  various 
levels  of  metal  content,  changes  in  viscosity,  and  failure 
of  a  blotter  test.  Baseline  data  provided  by  the  personnel 
at  the  JOAP  Laboratory  for  determining  failure  of  each  type 
of  oil  is  in  Appendix  0. 

The  entire  set  of  raw  data  was  provided  by  the  JOAP 
Laboratory  for  this  research.  To  maintain  an  unbiased 
sample,  we  collected  only  one  instance  of  oil  failure  per 
vehicle  in  the  test,  although  many  of  the  vehicles 
experienced  more  than  one  oil  change  recommendation.  To 
have  collected  more  than  one  lifetime  per  vehicle  would  have 
biased  our  analysis  in  favor  of  those  vehicles  with  short 
oil  lifetimes.  Data  was  compiled  in  terms  of  several 
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■factors  as  shown  at  Appendix  E. 

Although  we  received  the  raw  data  for  approximately 
450  Air  Force  vehicles,  a  thorough  examination  of  the  data 
revealed  several  problems  which  significantly  reduced  the 
total  number  in  our  research  sample.  For  example,  in  many 
cases,  an  accurate  number  of  miles  between  oil  changes  could 
not  be  determined.  Some  oil  was  not  changed  when  failure 
was  apparent  and  a  recommendation  was  made.  In  some  cases, 
the  wrong  type  of  oil  was  added  to  the  engine  crankcase, 
invalidating  the  test.  Also,  two  bases  were  removed  from 
the  project  due  to  engine  problems  and  the  use  of  improper 
oil.  After  discarding  the  samples  we  felt  were  inaccurate, 
the  total  number  of  cases  remaining  was  116.  While  this 
number  is  certainly  not  as  large  as  originally  anticipated, 
we  consider  it  a  fair  cross-section  of  oil  types  and  bases. 

Statistical  Analysis 

Prior  to  performing  the  TTT  analysis,  we  ran  several 
Ui 1 coxon-Rank  Sum  <Mann-Wh i tney )  tests  on  the  data  set. 

This  test  is  a  nonpar ame tr i c  technique  for  testing  the  null 
hypothesis  that  the  probability  distributions  associated 
with  two  populations  are  identical.  We  made  comparisons  for 
three  factorsi  oil  type,  vehicle  make,  and  utilization  and 
mileage.  Based  on  an  alpha  value  of  .05,  in  only  two  cases 
could  we  conclude  that  the  populations  were  significantly 
different.  The  failure  distribution  for  vehicles  having  low 
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utilization  and  low  mileage  was  significantly  different  from 
both  the  distribution  for  vehicles  having  high  utilization 
and  high  mileage  and  the  distribution  for  vehicles  having 
low  utilization  and  high  mileage.  These  two  tests  resulted 
in  two-tailed  probability  values  of  .’0097  and  .0105, 
respecti vely.  However,  there  was  no  significant  difference 
between  the  failure  distributions  based  on  oil  type,  vehicle 
make,  or  the  other  utilization  and  mileage  categories. 

Based  on  these  results,  we  eliminated  the  low 
utilization/low  mileage  vehicles  from  the  data  set  and 
retested  for  significant  differences  between  oil  types  and 
vehicle  makes.  Again,  we  were  unable  to  find  significant 
differences  among  any  of  the  failure  distributions. 

Finally,  we  tested  only  the  vehicles  having  low 
utilization/low  mileage  for  significant  differences  based  on 
oil  type  and  vehicle  make.  Again,  we  found  no  significant 
differences.  Thus,  based  on  our  sample  data,  we  have  no 
evidence  that  either  oil  type  or  vehicle  make  significantly 
affects  the  propensity  of  motor  oil  to  fail.  Results  of 
these  tests  are  shown  at  Appendix  F. 

Total  Time  on  Test 

To  demonstrate  the  Blind  Replacement  model  using 
graphical  analysis,  a  value  for  replacement  time  K  must  be 
established.  We  estimated  one  day  as  the  approximate  time 
required  to  turn  a  vehicle  into  the  motor  pool,  change  the 
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oil  and  filter,  and  return  it  to  service.  Since  the  life 
cycle  length  in  the  blind  replacement  model  is  expressed  as 
replacement  interval  plus  replacement  time,  (T+K) ,  T  and  K 
must  have  the  same  unit  of  measure.  Also,  since  it  is 
desirable  to  express  the  replacement  interval  in  miles, 
which  is  the  common  unit  of  measure  for  oil  life,  in  order 
to  maintain  consistent  units  of  measure,  K  should  also  be 
expressed  in  miles.  Therefore,  one  day  is  translated  into 
miles  by  dividing  the  total  miles  driven  in  the  test  by  the 
number  of  days  in  the  test  to  obtain  an  estimate  of  miles 
per  day.  Thus,  replacement  time  is  being  measured  as  the 
expected  number  of  miles  that  cannot  be  driven  when  the  oil 
is  being  replaced.  Using  this  technique,  the  following 
values  for  K  were  obtained: 

All  Oils  Combined  —  24.83 

Stauffer  Synthetic  Oil  *  20.92 
Mineral  Oi 1  »  28.30 

CONOCO  Synthetic  Oil  =  25.46 

Ue  developed  the  FORTRAN  computer  program  shown  at 
Appendix  0  to  perform  the  TTT  analysis.  Using  the  computed 
values  for  K,  we  applied  our  data  set  to  this  program  and 
found  the  following  replacement  intervals  provide  maximum 
availability  of  good  oil: 

All  oils  combined  (Appendix  H)  =»  1971  miles 
Stauffer  (Appendix  I)  *  1139  miles 

Mineral  oil  (Appendix  J)  *  2039  miles 

CONOCO  (Appendix  K)  =  2624  miles 
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Unfor  tunatel  y ,  in  all  four  cases,  the  value  of  -K/ t  ^  n  ^  jS  so 
small,  that  when  graphed,  it  caX  hard)/  be  distinguished 
•from  the  origin.  Thus,  the  TTT  plots  in  Appendices  H 
through  K  do  not  provide  very  good  illustrations  of  the  TTT 
process. 

Summary 

Based  on  the  tests  performed,  we  conclude  that  the 
failure  probability  distributions  are  not  significantly 
different  for  the  three  types  of  oils  tested.  However, 
since  the  TTT  transform  can  select  the  optimal  interval  from 
only  the  observed  data  points,  it  would  not  select  the  same 
replacement  interval  for  all  three  type  oils  unless  each 
happened  to  have  an  oil  change  point  of  the  same  mileage  and 
the  value  for  K  in  each  sample  allowed  that  interval  to  be 
selected.  For  this  reason,  any  interpretation  of  the 
differences  in  the  intervals  should  be  approached  with 
caution.  It  should  also  be  noted  that  the  optimal  interval 
selected  is  for  the  prescribed  objective  only — specifically, 
to  maximize  the  availability  of  good  oil.  If  the  decision 
maker  has  other  objectives,  these  intervals  would  no  longer 
be  optimal  since  use  of  the  blind  replacement  model  would  be 
i nappropr i ate . 
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CHAPTER  V 


4 


MODEL  ANALYSI S 

While  analyzing  the  results  of  this  project,  we 
noted  major  difficulties  in  our  application  of  the  Talbott 
model  and  weaknesses  in  the  original  Radner-Jorgenson 
formulation  in  terms  of  its  real -world  applicability.  In 
this  chapter,  we  will  describe  the  weaknesses  of  the  Blind 
Replacement  model,  explain  how  the  Talbott  model  was 
misapplied  for  our  data  set  and  then  present  a  proposed 
methodology  for  overcoming  some  of  these  problems. 

Blind  Replacement  Model 

The  weaknesses  of  the  blind  replacement  model  can 
best  be  described  in  terms  of  Markov  states.  As  stated 
earlier,  type  two  uncertainty  exists  when  there  is  doubt 
regarding  the  present  condition  of  the  device.  Thus,  as  a 
device  operates  under  conditions  of  type  two  uncertainty, 
can  be  in  any  one  of  four  possible  states! 

1.  Good  time  -  That  time  during  which  the  device  is 
in  operation  in  an  unfailed  condition. 

2.  Bad  time  -  That  time  during  which  the  device  is 
in  operation  in  a  failed  condition. 

3.  Corrective  maintenance  -  That  time  during  which 
a  failed  device  is  being  replaced. 

4.  Preventive  maintenance  -  That  time  during  which 
an  unfailed  device  is  being  replaced. 
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In  terms  of  the  lift  cycle,  at  any  point  in  tint* 
during  the  replacement  interval  T,  th*  device  is  in  either 
state  on*  or  state  two;  throughout  the  replacement  time  K, 
the  device  is  in  either  state  three  or  state  four.  Since  a 
device  can  be  available  only  when  it  is  in  state  one,  when  T 
equals  T* ,  the  probabi 1 i ty  of  being  in  state  one  versus  any 
other  state  is  maximized.  Hence,  blind  replacement 
emphasizes  state  one  at  the  expense  of  all  the  other  states 
combined.  Therefore,  the  probability  of  being  in  states 
two,  three,  and  four  is  minimized  unconstrained  with  respect 
to  each  other  as  if  they  were  a  single  state.  Since  there 
cannot  be  a  direct  transition  from  state  one  to  state  three, 
only  the  conditions  for  states  two  or  four  can  take  a  device 
out  of  state  one.  Thus,  blind  replacement  assumes  an  equal 
preference  for  entering  states  two  and  four  and  further 
assumes  that  they  both  have  the  same  intrinsic  value  in  the 
mind  of  the  decision  maker.  In  other  words,  one  must  be 
able  to  say,  "I  would  like  the  device  to  be  available;  but 
if  it  must  stop  being  available,  it  makes  no  difference 
whether  this  is  because  it  is  staying  in  operation  in  a 
failed  condition  or  because  it  is  undergoing  preventive 
replacement.*  Thus,  where  most  models  assume  a  cost  of 
failure  that  is  greater  than  the  cost  of  replacement,  the 
Radner- Jorgenson  model  implicitly  assumes  these  costs  are 
equal  or  of  no  concern. 
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Radner  and  Jorgenson  specifically  developed  their 
model  for  devices  (such  as  munitions)  being  Kept  in  storage 
under  conditions  of  type  two  uncertainty  for  use  during  an 
emergency.  Here,  the  objective  is  to  maintain  the  device  in 
a  state  of  operational  readiness  (state  one)  as  much  of  the 
time  as  possible.  Hence,  states  one  and  two  constitute 
times  that  the  device  is  in  storage,  not  in  operation.  For 
this  scenario,  there  may  be  little  difference  between 
failure  and  replacement  costs.  When  the  emergency  hits,  the 
cost  of  not  having  the  device  operationally  ready  is  the 
cost  of  not  being  in  state  one.  This  cost  is  incurred 
regardless  of  which  of  the  other  three  states  the  device  is 
inj  hence,  it  is  not  a  failure  cost.  The  only  failure  cost 
is  the  cost  of  the  type  two  uncertainty,  which  is  inherent 
in  the  problem  to  begin  with.  When  the  device  is  in  state 
three  or  four,  the  device  is  known  to  be  not  available  to 
respond  to  the  emergency.  Thus,  the  cost  of  failure  is  the 
cost  of  not  having  Known  the  device  to  be  failed. 

Where  this  information  has  appreciable  value,  the 
blind  replacement  model  should  not  be  applied.  Such  a 
situation  exists  in  the  case  of  a  munition  device  being 
"stored*  on  another  piece  of  operational  equipment,  such  as 
an  aircraft  ejection  seat  cartridge.  Here,  the  cost  of 
failure  is  possible  loss  of  life  since,  Knowing  the  device 
to  be  faileo,  one  would  not  fly  the  aircraft.  Th?  situation 
also  exists  for  munitions  on  alert  if,  for  planning 
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purposes,  it  is  considered  advantageous  to  reduce  type  two 
uncertainty  at  the  expense  of  availability. 

Another  problem  with  application  o-f  the  blind 
replacement  model  is  most  evident  when  there  is  a  very  small 
value  for  K  relative  to  t(n).  Necessarily,  as  K 
approaches  0,  T  approaches  t^jj.  In  fact,  if  K  **  0,  any 
value  for  T,  where  0  i  T  <.  t^^  will  provide  100 V.  avail¬ 
ability.  Therefore,  with  a  very  small  value  for  K,  maximum 
avai lability  is  realized  when  T  *  *(!)•  *n  this  situation, 
all  but  one  of  the  devices  will  be  replaced  while  they  still 
have  life,  and  this  could  be  costly.  Thus,  there  is  another 
cost  involved:  the  cost  of  discarded  life. 

In  a  minimum  cost  model,  not  accounting  for 
discarded  life  biases  the  model  toward  an  earlier 
replacement  interval.  The  greater  the  difference  between 

failure  and  replacement  costs,  the  relatively  less  expensive 

/• 

will  be  the  discarded  life,  but  the  more  discarded  life  will 
be  incurred.  However,  if  failure  and  replacement  costs  are 
relatively  close,  discarded  life  becomes  relatively 
expensive  but  rarely  occurs.  Thus,  when  minimizing  costs, 
even  though  the  bias  exists,  the  very  nature  of  the  model 
will  constrain  the  significance  of  discarded  life  costs.  On 
the  other  hand,  when  the  objective  of  the  model  is  to 
maximize  availability,  the  model  will  discard  life 
unconstrained.  If  the  device  is  relatively  inexpensive  and 
K  is  relatively  large,  discarded  life  costs  will  be  minimal. 
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However,  if  the  r»v#r*»  is  true,  discarded  life  could  be 
quite  significant,  and  the  costs  should  not  be  ignored. 

In  conclusion,  the  Radner-Jorgenson  and  the  Talbott 
models  appear  to  be  valid  when  applied  strictly  to  the 
scenario  for  which  they  were  developed;  however,  they  appear 
to  be  weak  in  terms  of  real-world  applications.  The 
inherent  problem  with  blind  replacement  is  not  in  the 
Radner-Jorgenson  formulation  but,  rather,  in  the  objective 
function.  Any  model  that  maximizes  availability  without 
constraint  is  not  a  very  useful  model  in  our  opinions.  To 
the  extent  that  there  is  a  preference  for  entering  state 
four  versus  state  two  or  that  discarded  life  is  a  concern, 
maximizing  availability  without  constraint  is  not  realistic. 
In  other  words,  costs  must  truly  be  irrelevant  for  this 
objective  to  make  sense,  and  we  doubt  this  is  often  the 
case . 

Appl  iSiUQD  9*  J&UlgH 

The  application  of  the  Talbott  model  in  this 
research  project  for  finding  an  optimal  oil  replacement 
interval  was  inappropriate  since  engine  oil  does  not 
constitute  a  device  in  storage  waiting  to  respond  to  an 
emergency,  and  there  are  real  costs  associated  with  failure 
and  discarded  life  which  probably  should  not  be  ignored. 
Hence,  the  optimal  replacement  interval  is  only  valid  if  the 
objective  is  to  maximize  avai 1 i abi 1 i ty  of  the  good  oil  and 
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there  is  no  concern  -for  the  cost  of  operating  an  engine  on 
abada  oil  or  for  the  cost  of  discarding  good  oil  due  to 
early  replacement. 

Uhat  would  be  an  appropriate  model  for  oil 
replacement  is  a  question  the  authors  are  unable  to  answer. 
There  are  several  problems  associated  with  engine  oil  that 
make  application  of  an  appropriate  replacement  model 
complex.  First,  as  oil  transitions  from  state  one  to  state 
two,  failure  costs  are  not  incurred  all  at  once;  they  are 
accumulated  as  the  engine  continues  to  operate  on  failed 
oil.  Second,  this  accumulation  of  failure  costs  is  probably 
not  linear  since  operating  an  engine  on  failed  oil  will 
probably  aggravate  engine  damage  in  a  nonlinear  fashion. 
Thus,  the  cost  of  operating  a  vehicle  for  a  mile  on  oil  that 
Just  failed  is  probably  less  than  the  cost  of  operating  a 
vehicle  for  a  mile  on  oil  that  has  been  failed  for  a  long 
time.  Finally,  this  process  of  how  failed  oil  damages  an 
engine  over  time  is  probably  stochastic  and  is,  therefore, 
represented  by  some  theoretical  probability  distribution. 

So,  even  if  one  were  to  obtain  an  unbiased  sample  of  this 
process,  the  issue  of  how  to  address  the  type  three 
uncertainty  associated  with  this  theoretical  distribution 
would  still  be  a  problem.  Thus,  considering  the 
intractabi 1 i ty  of  accounting  for  failure  cost,  proposing  an 
appropriate  model  for  determining  an  optimal  oil  replacement 
interval  presents  a  dilemna  the  authors  have  been  unable  to 
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resol ve 


PrgB.gt.tl  OtlHadfll  asm 

As  stated  earlier,  type  three  uncertainty  exists  if 
the  underlying  failure  distribution  is  not  known  with 
certainty.  Both  models  that  have  been  developed  to  address 
this  type  uncertainty,  Bergman's  and  Talbott's,  are  somewhat 
limited  in  terms  of  real-world  applications  since  neither 
considers  the  cost  of  discarded  life  and  both  maximize  an 
objective  function  without  constraint.  The  Bergman  model 
minimizes  cost  without  regard  for  availability,  and  the 
Talbott  model  maximizes  availability  without  regard  for 
cost.  Therefore,  the  authors  propose  two  models  termed 
model  one  and  model  two  which  allow  the  decision  maker  to 
constrain  the  objective  function.  Both  repl icemen t  models 
one  and  two  address  three  concerns  which  the  authors 
consider  appropriate  for  most  real-world  scenarios.  These 
concerns  are  cost,  availability,  and  risk  of  failure.  When 
the  decision  maker  selects  one  of  these  concerns  as  the 
objective,  the  other  two  may  remain  in  the  model  as 
constraints.  We  will  address  replacement  model  two  first. 

Replacement  model  two  is  appropriate  for  situations 
having  all  three  types  of  uncertainty.  To  illustrate 
application  of  this  model,  we  recall  the  example  of  the 
aircraft  ejection  seat  cartridge.  An  appropriate 
application  in  this  example  might  be  to  minimize  cost, 
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constrained  by  an  acceptable  level  of  risk  -for  being  in 
operation  in  a  -failed  state  (state  two).  Furthermore,  in 
the  example  o-f  munitions  on  alert,  it  might  be  more 
realistic  to  address  availability  as  a  constraint,  not  an 
objective.  Thus,  an  appropriate  application  might  be  to 
minimize  cost,  constrained  by  a  required  level  o-f  readiness 
to  meet  mission  requ i remen ts .  Here,  where  availability 
operates  as  a  constraint,  the  cost  o-f  not  having  known  the 
device  to  be  failed  is  minimized  over  the  long  run.  By 
knowing  the  availability  constraint,  the  decision  maker 
always  knows  the  expected  number  of  devices  in  the  readiness 
state  (state  one)  and  can  plan  accordingly. 

In  the  model  two  formulation,  T  represents  the  re¬ 
placement  interval;  C(T)  is  the  expected  long  term  cost; 

Cj  is  the  cost  of  device  failure  to  include  cost  of 
replacing  the  failed  device;  and  Cg  is  the  cost  of  replac¬ 
ing  a  good  device.  Since  discarded  life  is  the  cumulative 
life  within  the  interval  (T,oo),  the  expectation  for  the 
percent  of  the  total  life  discarded  can  be  expressed  as: 

oo 

(  R(x)dx 

i 

E(x)  * -  (5.1) 

(°° 

)  R(x)dx 

0 

LJhere  cost  of  discarded  life  is  measured  as  the  cost  of  re¬ 
placement  multiplied  by  the  expected  percent  of  the  total 


life  discarded,  then  cost  of  discarded  life  equals 
oo  oo 

C?C  (  R<x)dx/  )  R(x)dx).  Further,  in  accordance  with 
TT  0 

renewal  reward  theorem,  the  objective  of  a  minimum  cost 
model  is  to  minimize  expected  costs  per  cycle  divided  by  the 
expected  cycle  length  (i.e.,  minimize  cost  per  unit  of  life 
such  as  dollars  per  operating  hour).  Hence,  similar  to 
Barlow  and  Proschan's  age  replacement  model,  the  optimum 
replacement  interval  for  minimizing  cost  is  that  value  of  T 
which  yields  the  minimum  value  for  C(T)  in  the  following 
equat i on : 


C  F<T)  ♦  C0R<TH  \  R(x)dx/  \  R(x)dxl 
1  2  t'  O' 

C(T)  * -  (5.2) 

T 

<  R(x)dx 

o' 


For  maximizing  avai 1 abi 1 i ty,  the  optimum  replacement 
interval  is  found  using  the  Radner-Jorgenson  formulation, 
which  maximizes  the  expected  life  over  the  expected  cycle 
time  of  T  ♦  Kj  specifically, 


T 

<  R(x)dx 

o' 

A<T)  - -  (5.3) 

T  ♦  K 


For  minimizing  risk  of  failure,  we  let  Z(T) 
represent  expected  long  term  risk  of  the  device  being  in 
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state  two.  Since  the  objective  is  to  minimize  the  expected 
time  in  a  -failed  condition  over  the  expected  cycle  time,  the 
optimum  replacement  interval  is  that  value  o-f  T  which 
minimizes  Z(T>  in  the  -following  equation, 

T 

(  F<x)dx 

0 

Z  (T)  * -  (5.4) 

T  +  K 

Minimizing  Z(T)  unconstrained  would  not  make  sense  since  it 
would  necessarily  drive  a  replacement  interval  o-f  T  —  0, 
where  the  expected  risk  would  also  be  0. 

In  order  for  the  model  to  address  type  three 
uncertainty  (i.e.,  uncertainty  regarding  the  underlying 
failure  distribution),  we  apply  some  lessons  learned  from 
total  time  on  test  methodology.  As  stated  earlier, 

represents  a  data  point  in  the  ordered  set  of  life 
observations;  T.  represents  the  accumulated  life  at 
point  in  time  t<f);  and  n  represents  the  number  of 
observations.  Then,  we  let  £.  represent  a  cost  per  unit  of 
life  constraint;  z  represent  an  acceptable  level  of  risk  for 
the  device  being  in  state  two;  a  represent  an  availability 
constraint  in  terms  of  the  probability  of  the  device  being 
in  state  one;  B.  represent  the  accumulated  time  in  state 
two;  and  0.  represent  the  accumulated  life  discarded  by 
preventive  replacement.  Then  for  any  value  of  T  =  *<j)> 
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i-1 


j-1 


(5.5) 


n 

°i  -  H  <*<j>  •  *<(>’ 

j«i  +  l 


<5. 6) 


T+K  -  t(  j+K 


<5.7> 


Since  T./n  equals  the  accumulated  life  for  the  n  components 
divided  by  the  number  of  components,  it  necessarily  equals 
the  expected  life  for  the  empirical  data  set  and  will 
approximate  the  mean  of  the  theoretical  probability 
distribution  when  n  is  sufficiently  large.  Thus,  for  any 
val ue  of  T  »  t( . ^  , 

3  1 im  T./n  <5. 8) 

n — >oo  i 

-nJ4&Vn  <5‘9> 


T 

$  R( x ) dx 


And  si  mi  1 arl y , 


.00 


j  R<x)dx 


oo 

R<x)dx  —  1 im  D./n 

n— >oo  i 


(  F  <  x  >  dx 

0 


1 im  B./n 
n— >oo  i 


<3. 10) 

(5.11) 
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Also,  since  F(T)  =  1 im  F _<T) ,  then  when  n  is  sufficiently 

n— >oo 

1  arge ,  F<T>  4  Fn<T>  -  Fn(t<i>)m  And 

R<T>  4  Rn(t(|))|  A<T>  *  An<t(j>)j  C<T>  4  Cn<t(j)>l  *nd 
2<T>  4  Zn<t<j>>*  Further,  since  t(j)  represents  one  data 

point  in  the  empirical  cumulative  failure  distribution, 
Fn<T>,  then  necessarily  Fn<t<j)>  i*  the  probability  of 
failure  within  the  interval  *t<  j  *n<*  's  therefore 

approx imatel y  equal  to  i/n.  Similarly,  Rn<t^jj>  4  l-<i/n>. 
Thus  from  equation  5.2, 


C. < i/n)  +  C~Cl-< i/n)]<D;/T  > 
C  <t/iN)  *  — 1 - 2 - l—Q 


n  <i> 


T./n 


(5.12) 


C^< i )  +  Cgtn-i ) <Dj/Tn) 


And  from  equation  5.3, 


T./n 

A  (tn  J  ^  - L - 

n  ( i  >  *  +  k 

t(i>  +  K 


-li. 


n(t( j  j+K) 


(5.13) 


And  from  equation  5.4, 


v  .  B;/n 

Z  <  t ,  .  .  >  =  - + - 

n  ( •  >  *■  a  i/ 

t(j)  ♦  K 


_ 

n(t( . j+K) 


(5.14) 


Therefore,  we  can  conclude  that! 
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min  C<T>  “> 


max  A<T> 


min  Z<T) 


min 


Cjd)  ♦  C2<n-i  ><Dj/Tn> 
T. 


>  max 


n<t( , }+K> 


B. 


«>  min 


n<t( . j+K) 


<3.13) 


<5. 16) 


<5.17) 


Thus,  the  value  of  t(()  that  minimizes  C<T>  where 
Bj/[n<t< . j+K))  <  z  and  Tj/En< t( f j+K) ]  >  a  is  tha  replacement 
interval  that  minimizes  long  term  costs  subject  to  the 
probability  o-f  the  device  being  in  state  two  not  exceeding  z. 
and  probability  o-f  tha  device  being  in  state  one  being 
greater  than  or  equal  to  j,.  Likewise  the  value  o-f  t(i)  that 
maximizes  A<T)  where  Zn<t(j))  <.  j,  and  Cn<t(j))  <.  c  i  s  the 
replacement  interval  that  maximizes  availability  subject  to 
probability  of  the  device  being  in  state  two  not  exceeding  z. 
and  cost  per  unit  of  life  not  exceeding  Either  of  these 
solutions  can  be  found  using  a  simple  iterative  computer 
program.  Also,  note  that  finding  the  value  of  t^ . ^  that 
maximizes  A<T)  where  there  are  no  constraints  established 
for  cost  and  risk  of  failure  is  equivalant  to  finding  the 
optimal  replacement  interval  using  the  Talbott  model.  Thus, 
we  now  formulate  replacement  model  one. 

Replacement  model  one  is  appropriate  for  situations 
having  only  type  one  and  type  three  uncertainty.  The  major 
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distinction  in  this  scenario  is  that  since  there  is  no  type 
two  uncertainty  <i.e.t  uncertainty  regarding  the  present 
condition  of  the  device),  a  device  would  not  stay  in 
operation  in  a  -failed  condition;  there-fore,  there  is  no 
state  two.  Hence,  the  optimum  replacement  inteval  -for 
maximizing  availability  is  that  value  o-f  T  which  maximizes 
the  expected  life  over  the  expected  cycle  time  in  the 
following  equation: 


T 

(  R<x)dx 

o' 

ACT)  =  -  (5.18) 

T 

(  R<x)dx  +  K 

0 


Here,  maximizing  ACT)  unconstrained  would  not  make  sense 
since  it  would  necessarily  drive  a  replacement  interval  of  T 
equal  to  infinity.  For  minimizing  risk  of  failure,  the 
optimum  replacement  interval  is  simply  that  value  of  T  that 
minimizes  FCT) .  The  minimum  cost  formulation  is  identical 
to  the  formulation  for  replacement  model  two  CEquation  5.2). 
Therefore,  we  can  conclude: 


min  C<T> 


min 


Cl<i>  +  CgCn-i )COj/Tn) 


C5.19) 


max  ACT) 


max 


T. 

i 


T  j  ♦  nK 


C5.20) 


48 


min  F<T>  ■■>  min  i/n  <5.21) 

Thus,  the  value  of  t(  js  that  minimizes  C(T>  where 
Tj/CTj+nK)  2  *  and  F<T>  i  j  is  the  replacement  interval  that 
minimizes  long  term  costs  subject  to  the  probability  o-f  the 
device  being  in  state  one  not  being  less  than  a  and  the  risk 
o-f  failure  not  exceeding  £.  Similarly,  maximum  availability 
and  minimum  risk  solutions  can  be  found  as  well.  Further¬ 
more,  finding  the  value  of  t^fj  that  minimizes  C(T>  where 
i  »  1  and  £  »  0  is  similar  to  finding  the  optimal 
replacement  interval  using  the  Bergman  model.  The  only 
d i-fference  is  that  the  Bergman  model  does  not  account  for 
discarded  life. 

Ssfflnarz 

Achieving  the  best  value  for  a  prescribed  objective 
function,  unconstrained  by  any  other  factors,  may  not  be  a 
realistic  goal  for  many  real-world  preventive  maintenance 
scenarios.  This  is  particularly  true  when  availability 
(readiness)  or  risk  of  failure  are  viable  concerns,  as  they 
often  are  in  the  Air  Force.  Under  conditions  of  type  two 
uncertainty  <i.e.,  uncertainty  regarding  the  present 
condition  of  the  device),  maximizing  availability  without 
constraint  can  lead  to  several  problems,  as  discussed 
earlier)  and  minimizing  risk  of  failure  without  constraint 
may  not  be  prudent.  In  the  absence  of  type  two  uncertainty, 
neither  maximizing  availability  nor  minimizing  risk  of 
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■failure  makes  any  sense  as  an  unconstrained  objective 
-function.  Hence,  we  believe  models  that  allow  the  decision 
maker  to  impose  constraints  on  the  objective  function  are 
far  more  realistic. 

Identifying  an  appropriate  model  for  selecting  an 
optimal  oil  replacement  interval  remains  a  dilemna.  Due  to 
the  complex  nature  of  motor  oil  failures,  none  of  the  models 
addressed  in  this  project  appear  to  be  appropriate.  We 
believe  selecting  an  optimal  preventive  maintenance  interval 
for  motor  oil  requires  development  and  application  of  a 
sophisticated  maintenance  model,  far  more  complex  than  the 
ones  addressed  in  this  project. 
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CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  purpose  of  this  chapter  is  to  draw  conclusions 
from  this  research  project  and  make  some  recommendations  for 
further  research  and  study. 

Research  Question  and  Conclusions 

The  research  question  for  this  project  was,  "Can  the 
blind  replacement  model  be  applied  to  an  empirical  data  set 
using  TTT  transform  graphical  solution  techiques  to 
determine  optimal  replacement  intervals  for  synthetic  and 
petroleum  motor  oils?"  Our  answer  is  no.  The  major  reason 
is  that  motor  oils  do  not  meet  the  assumptions  of  the  blind 
replacement  model  since  they  are  not  in  storage  waiting  to 
respond  to  an  emergency.  Also,  we  believe  there  are  viable 
cost  concerns  associated  with  replacing  motor  oils. 
Specifically,  there  are  real  costs  resulting  from  operating 
an  engine  on  failed  oil  and  discarding  life  of  good  oil. 
Moreover,  these  costs  probably  should  not  be  ignored  when 
selecting  oil  replacement  intervals.  Since  blind 
replacement  maximizes  availability  without  regard  for  cost, 
we  consider  it  an  inappropriate  model  for  determining 
optimal  replacement  intervals  for  motor  oils.  Therefore,  we 
cannot  recommend  the  replacement  intervals  determined  in 
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this  research  project  for  real-world  applications. 

In  addition  to  the  above,  the  authors  draw  the 
-following  conclusions  based  on  the  data  and  analysis  of  this 
research  project. 

1.  When  comparing  synthetic  (CONOCO  and  Stauffer) 
and  petroleum  (mineral)  motor  oils,  the  type  of  oil  used  in 
an  Air  Force  general  purpose  vehicles  does  not  significantly 
affect  the  propensity  of  the  oil  to  fail.  In  other  words, 
there  is  no  significant  difference  between  these  three  types 
of  motor  oils  in  terms  of  probability  of  oil  failure. 

2.  In  our  opinions,  due  to  the  complex  nature  of 
oil  failures,  application  of  an  appropriate  model  for 
determining  optimal  replacement  intervals  requires  selection 
of  a  fairly  complex  preventive  maintenance  model. 

3.  Both  the  Radnei — Jorgenson  blind  replacement 
model  and  the  Talbott  adaptation  of  this  model  to  TTT 
transform  techniques  appear  to  be  valid  when  applied 
strictly  within  the  scenario  for  which  they  were  developed. 
However,  we  believe  their  objective  of  maximizing 
availability  without  constraint  makes  them  inappropriate  for 
many  real-world  applications. 

4.  The  TTT  transform  technique  appears  to  be  a 
viable  approach  for  determining  optimal  replacement 
intervals  when  the  underlying  failure  distribution  is  not 
known  with  certainty  since  it  eliminates  error  resulting 
from  type  three  uncertainty. 
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5.  Preventive  maintenance  models  that  achieve  the 
best  value  of  a  prescribed  objective  function,  unconstrained 
by  any  other  factors,  may  not  be  appropriate  for  many 
real-world  preventive  maintenance  scenarios. 

Recommendations  for  Further  Research 

The  authors  recommend  the  following  areas  for  future 
research  and  study: 

1.  Development  and  demonstration  of  preventive 
maintenance  models  for  motor  oils  and  devices  having  similar 
complex  failure  patterns. 

2.  Further  adaptations  of  mathematical  models  to 
TTT  transform  techniques  or  other  methodologies  that 
eliminate  error  resulting  from  type  three  uncertainty. 

3.  Further  development  and  demonstration  of 
preventive  maintenance  models  that  allow  constraints  to  be 
imposed  on  the  objective  function. 

4.  Greater  focus  in  the  arena  of  reliability  theory 
toward  demonstration  of  preventive  maintenance  models  in 
terms  of  their  real-world  tractability  and  applications. 
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TTT  PROCESS  FOR  AN  IFR  DISTRIBUTION 


EXAMPLE :  Consider  a  set  of  five  components  where  one  fails 

after  3  hours  of  operation,  one  after  9  hours,  one  after  12 
hours,  one  after  14  hours,  and  one  after  13  hours.  The 
ordered  set  of  life  distributions  is  then  expressed  as 
(3,  9,  12,  14,  15).  From  equation  3.1,  T.  is  computed  for 
each  component  where  n  »  5  as  follows: 
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T. 
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-n 

<5  - 

■  j 

♦ 

1)  <t(j)  -  t(j_ 

1)> 

j*l 
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*  54 

T5 
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♦ 

1) 

<15  -  14)  ♦  54 

*  55 

And  Uj  can  be  computed  as  Tj/Tn;  hence, 
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at  Attachment  C. 
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APPENDIX  B 

EXiVIPLE  TTT  PROCESSES  FOR  EXPONENTIAL  AND  DFR  DISTRIBUTIONS 
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EXAMPLE  TTT  PROCESSES  FOR  EXPONENTIAL  AND  DFR  DISTRIBUTIONS 


Example  Exponential  distribution! 
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Ex  amp  I  e 

DFR  Distribution: 
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TTT  plots  for  these  distributions  are  given  at  Attachment  C 


APPENDIX  C 
EXAMPLE  TTT  PLOTS 
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EXAMPLE  TTT  PLOTS 


Exponential  Distr 
(From  Appendix  B) 
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APPENDIX  D 

CRITERIA  FOR  FAILED  OIL 
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SPtCTROMtTR IC  ANALYSIS 
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6.02 

M 

Total  Base  Number-D2S96 

mgKOH/ g 

6.48 

6.32 

It 

Total  Acid  Number- D664  , 

mgKOH/g 

2.75 

2.81 

•  < 

Sulfated  Ash,  l 

1.02 

0.99 

0.90-1.20 

CRC  L-38 

Crankcase  Oil  Oxidation 

Test 

Bearing  weight  loss  (mg) 

— 

Top  half 

13.5 

16.2 

- 

Bottom  half 

17.4 

15.9 

- 

Total 

30.9 

32.1 

<40.0 

Varnish  Deposits 

59.3 

59.4 

- 

Piston  Skirt  Varnish 

9.6 

9.6 

V  9 . 0 
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CONOCO  OM-oOO 


3126 


GENERAL  PURPOSE  t$~xjo* i 
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Vehicle  make: 

AMC  =>  American  Motors  Corporation 
CHEV  -  Chevrolet 
0006  *  Dodge 
FORD  «  Ford 

IHC  »  International  Harvester  Corporation 
OLDS  91  Oldsmobile 
PLYM  *  Plymouth 

Vehicle  Type: 

AMB  =  Ambulance 

C/A  ■  Carry-all 

CJ7  =»  Jeep 

M/S  =  Multi-step 

P/U  *  Pick-up 

Pan  *  Panel  truck 

S/P  m  Special  purpose 

S/V  *  Step  van 

SDN  =  Sedan 

SU  =  Station  wagon 

TRAC  “  Tractor 

1  Ton  =  One  ton  truck 

1.3TN  «  1-1/2  ton  truck 

4X4  »  Four-wheel  drive  truck 

Usage  Code : 

A  ■  High  utilizati on/1 ow  mi  1  cage 
B  ”  High  utilization/high  mileage 
C  *  Low  utilization/low  mileage 
D  »  Low  utilization/high  mileage 
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LUBRICANT  LIFETIME  DATA 


STAUFFER 


MAKE  MODEL  USE  VEHICLE  ENG 


GRAND  FORKS  AF8 


CHEV 

C/A 

B 

791737 

350 

OFFUTT  AFB 

MINOT  AFB 

LACKLAND 

AFB 

one 

SDN 

D 

762826 

232 

DODG 

P/U 

C 

764048 

225 

fine 

S U 

A 

785344 

200 

fine 

SDN 

D 

788302 

232 

RANDOLPH 

AFB 

IHC  TRAC 

B 

724309 

392 

IHC 

M/S 

A 

732296 

345 

FORD 

P/U 

A 

761781 

1900 

AMC 

SU 

B 

762448 

232 

DODG 

P/U 

C 

770 1 49 

318 

CHEV 

P/U 

B 

781244 

250 

CHEV 

P/U 

C 

781247 

250 

FORD 

SU 

B 

785346 

200 

AMC 

SDN 

B 

788305 

232 

PLYM 

SON 

B 

789172 

225 

GEORGE  AFB 

CHEV 

S/V 

C 

761615 

350 

IHC  1 

.  5TN 

A 

771097 

345 

DODG 

P/U 

A 

780446 

318 

DODG 

P/U 

A 

780452 

318 

CHEV 

PAN 

A 

783521 

250 

FORD  1 

.  5TN 

A 

785850 

300 

PLYM 

SDN 

B 

789080 

225 

DODG  1 

TON 

fit 

789357 

360 

DODG 

FW 

A 

792534 

318 

START 

MILES 

END 

MILES 

OIL 

LIFE 

TEST 

DAYS 

44516 

48512 

3996 

85 

28114 

34081 

5967 

220 

23286 

30259 

6973 

402 

50426 

56651 

6225 

193 

10535 

25535 

15000 

442 

83212 

89272 

6060 

310 

54835 

60467 

5632 

204 

25220 

33194 

7974 

422 

39469 

48384 

8915 

305 

4544 

5251 

707 

330 

39335 

42231 

2896 

294 

13779 

19063 

5284 

440 

39725 

48238 

8513 

224 

16493 

26925 

10432 

361 

48537 

54752 

6215 

84 

27175 

32709 

5534 

339 

10003 

15554 

5551 

244 

23216 

26309 

3093 

187 

14853 

18137 

3284 

245 

29246 

34569 

5323 

289 

17445 

21467 

4022 

200 

33106 

38055 

4949 

93 

11639 

18660 

7021 

336 

11396 

13367 

1971 

134 

AVG 

MPD 


47.01 


27.12 

17.35 

32.25 

33.94 


19.55 

27.61 

18.90 
29.23 

2.14 

9.85 

12.01 

38.00 

28.90 
73.99 


16.32 

22.75 

16.54 

13.40 

18.42 

20.11 

53.22 

20.90 

14.71 
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MYRTLE  BEACH  AFB 


fiMC 

AMB 

C 

760278 

360 

2445 

3584 

1139 

527 

2.16 

DODG 

PAN 

A 

792473 

225 

2604 

11413 

8809 

586 

15.03 

HANCOCK 

fine 

FIELD 

SDN 

A 

733165 

304 

55928 

59672 

3744 

375 

9.98 

CHEV 

P/U 

C 

740738 

350 

58785 

62747 

3962 

385 

10.29 

FORD 

P/U 

A 

751429 

300 

35068 

43936 

8868 

497 

17.84 

FORD 

P/U 

A 

751430 

300 

29419 

33403 

3984 

381 

10.46 

PETERSON 

CHEV 

FIELD 

P/U 

A 

742105 

350 

46519 

52898 

6379 

320 

19.93 

DODG 

P/U 

B 

780139 

225 

23147 

28117 

4970 

209 

23.78 

DODG 

P/U 

C 

780141 

225 

7639 

11715 

4076 

383 

10.64 

CHEV 

M/S 

B 

784570 

292 

10645 

13244 

2599 

267 

9.73 

PLYM 

SW 

B 

785110 

225 

49242 

56395 

7153 

62 

115.37 

USAF  ACADEMY 
PLYM  SDN 

A 

788159 

225 

15345 

28276 

12931 

233 

55.50 

FORD 

SDN 

B 

795660 

200 

12535 

29364 

16829 

241 

69.83 
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MINERAL 


MAKE 

MODEL 

USE 

VEHICLE  ENG 

START 

MILES 

END 

MILES 

OIL 

LIFE 

TEST 

DAYS 

AVG 

MPD 

GRAND  FORKS  AFB 
CHEW  C/A  L 

791753 

350 

43211 

49814 

6603 

62 

106.50 

OFFUTT  AFB 
CHEV  M/S 

A 

722544 

350 

42094 

44133 

2039 

128 

15.93 

CHEV 

M/S 

A 

753604 

350 

41759 

49110 

7351 

195 

37.70 

MINOT  AFB 

CHEW  C/A 

C 

740088 

454 

164995 

170078 

5083 

403 

12.61 

CHEV 

C/A 

C 

744124 

454 

163183 

148057 

4874 

252 

19.34 

CHEV 

C/A 

B 

760519 

350 

1 17547 

133079 

15512 

214 

72.49 

CHEV 

C/A 

B 

791755 

350 

37750 

43697 

5947 

86 

69.15 

CHEW 

C/A 

B 

791764 

350 

49380 

53237 

3857 

20 

192.85 

CHEV 

C/A 

B 

792709 

350 

16038 

24616 

8578 

184 

46.62 

CHEW 

C/A 

B 

793380 

350 

17765 

28147 

10382 

149 

69.68 

CHEV 

C/A 

B 

793381 

350 

27407 

39101 

11694 

158 

74.01 

CHEV 

C/A 

B 

793383 

350 

30919 

39775 

8856 

116 

76.34 

LACKLAND 

DODG 

AFB 

C/A 

C 

755338 

318 

64995 

68712 

3717 

126 

29.50 

CHEV 

AMB 

A 

760113 

350 

49704 

55002 

5298 

405 

13.08 

AMC 

SDN 

B 

762919 

232 

45238 

53343 

8105 

267 

30.36 

CHEV 

P/U 

A 

781259 

250 

22498 

30926 

8428 

269 

31  .33 

DODG 

P/U 

C 

782261 

225 

17169 

22578 

5409 

247 

21  .90 

AMC 

SDN 

C 

788418 

232 

5603 

11205 

5602 

470 

11.92 

CHEV 

P/U 

B 

794279 

2300 

9171 

17362 

8191 

298 

27.49 

RANDOLPH 

CHEV 

AFB 

M/S 

A 

734641 

350 

52768 

56410 

3642 

140 

26.01 

CHEV 

C/A 

B 

770785 

350 

43730 

47060 

3330 

113 

29.47 

CHEV 

P/U 

B 

781245 

250 

19425 

35145 

15720 

506 

31  .07 

CHEV 

P/U 

B 

781246 

250 

28093 

35474 

7381 

371 

19.89 

FORD 

SW 

B 

787738 

200 

37805 

43587 

5782 

89 

64.97 

AMC 

SDN 

A 

788657 

232 

15792 

28178 

12386 

431 

28.74 
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GEORGE  AFB 


CHEV 

S/V 

B 

761619 

350 

36339 

39661 

3322 

170 

19.54 

CHEV 

S/V 

B- 

761622 

350 

41330 

46108 

4778 

147 

32.50 

IHC 

1  .  5TN 

A 

771098 

345 

29456 

35872 

6416 

206 

31 .15 

DODG 

P/U 

A 

780447 

318 

23776 

29663 

5887 

222 

26.52 

FLYM 

SDN 

B 

789081 

225 

34206 

37540 

3334 

133 

25.07 

DODG 

1  TON 

A 

789358 

318 

7151 

14730 

7579 

337 

22.49 

MYRTLE  BEACH  AFB 
OLDS  AMB  C 

730453 

455 

59300 

66539 

7239 

258 

28.06 

CHEV 

M/S 

A 

742721 

250 

41526 

46168 

4642 

511 

9.08 

DODG 

P/U 

C 

754213 

318 

14264 

17984 

3720 

276 

13.48 

HANCOCK  FIELD 


PETERSON  FIELD 


CHEV 

M/S 

B 

722632 

350 

53882 

59038 

5156 

235 

21  .94 

FORD 

P/U 

C 

761807 

1900 

31590 

41858 

10268 

387 

26.53 

DOCS 

4X4 

B 

770077 

225 

4708 

10049 

5341 

357 

14.96 

PLYM 

SU 

B 

785109 

225 

40787 

55673 

14886 

558 

26.68 

FORD 

SDN 

B 

783768 

400 

40832 

46272 

5440 

91 

59.78 

USAF  ACADEMY 

CHEV 

S/P 

A 

683520 

250 

62503 

66641 

4138 

299 

13.84 

DODG 

P/U 

A 

780131 

225 

19908 

25351 

5443 

240 

22.68 

DODG 

P/U 

A 

790861 

225 

6773 

18544 

11771 

231 

50.96 

CONOCO 


MAKE 

MODEL 

USE 

VEHICLE 

ENG 

START 

MILES 

END 

MILES 

OIL 

LIFE 

TEST 

DAYS 

AVG 

MPD 

GRAND  FORKS  AFB 

OFFUTT  AFB 

MINOT  AFB 

CHEV  C/A  C 

744121 

454 

159156 

164482 

5326 

293 

18.18 

CHEV 

C/A 

B 

760516 

350 

133095 

139086 

5991 

106 

56.52 

CHEV 

C/A 

B 

791757 

350 

56977 

61031 

4054 

54 

75.07 

CHEV 

C/A 

B 

791762 

350 

36924 

40912 

3988 

82 

48.63 

CHEV 

C/A 

B 

793387 

350 

14252 

20337 

6085 

308 

19.76 

CHEV 

C/A 

B 

793388 

350 

31466 

37194 

5728 

88 

65.09 

LACKLAND 

DODG 

AFB 

C/A 

A 

750870 

318 

51827 

69226 

17399 

602 

28.90 

DODG 

P/U 

D 

764049 

225 

37162 

47983 

10821 

267 

40.53 

CHEV 

P/U 

A 

781235 

250 

28602 

38397 

9795 

351 

27.91 

CHEV 

P/U 

A 

781261 

250 

12671 

23916 

11245 

428 

26.27 

AMC 

SU 

B 

785432 

200 

56456 

65300 

8844 

138 

64.09 

RANDOLPH 

FORD 

AFB 

P/U 

B 

761774 

1900 

46601 

56756 

10155 

216 

47.01 

CHEV 

C/A 

B 

770428 

350 

65402 

70377 

4975 

116 

42.89 

CHEV 

P/U 

B 

781248 

250 

17146 

25590 

8444 

318 

26.55 

FORD 

SU 

B 

785347 

200 

52692 

61004 

8312 

270 

30.79 

PLYM 

SDN 

B 

789173 

225 

63092 

75892 

12800 

147 

87.07 

GEORGE  AFB 
DODG  P/U 

B 

780442 

318 

30767 

33391 

2624 

188 

13.96 

DODG 

P/U 

B 

780460 

318 

27207 

30688 

3481 

122 

28.53 

FORD 

SU 

B 

787944 

302 

26953 

32057 

5104 

164 

31 .12 

PLYM 

SDN 

B 

789079 

225 

33043 

39500 

6457 

139 

46.45 

MYRTLE  BEACH  AFB 
DODG  P/U  C 

754212 

318 

14202 

181  18 

3916 

472 

8.30 

AMC 

CJ7 

C 

789616 

258 

16862 

21165 

4303 

322 

13.36 

W*JCOCK 

FORD 

FIELD 

SDN 

B 

714579 

302 

94743 

98742 

3999 

112 

35.71 

AMC 

SDN 

A 

733260 

304 

49598 

57033 

7435 

633 

11.75 

FORD 

P/U 

A 

751427 

300 

33875 

37462 

3587 

220 

16.30 

DODG 

P/U 

A 

754106 

318 

26777 

33360 

6583 

577 

11.41 

CHEV 

P/U 

A 

771056 

300 

12532 

18636 

6104 

545 

11.20 

PLYM 

SU 

B 

785038 

225 

19386 

27894 

8508 

152 

55.97 
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PETERSON 

CHEV 

FIELD 

M/S 

B 

722548 

350 

AMC 

SDN 

B 

733766 

304 

FORD 

P/U 

A 

753204 

300 

CHEV 

M/S 

B 

784568 

292 

CHEV 

M/S 

A 

784569 

292 

USAF  ACADEMY 
IHC  M/S 

A 

731969 

345 

FORD 

S/P 

A 

750245 

318 

DODG 

P/U 

A 

790858 

225 

FORD 

SDN 

B 

795659 

200 

56406 

63370 

6964 

296 

23.53 

47214 

55909 

8695 

495 

17.57 

39135 

44267 

5132 

552 

9.30 

19114 

27352 

8238 

278 

29.63 

18250 

24779 

6529 

231 

28.26 

64510 

71333 

6823 

323 

21.12 

43688 

47660 

3972 

184 

21  .59 

6934 

14413 

7479 

177 

42.25 

12161 

18434 

6273 

95 

66.03 

TOTAL  OIL  LIFE:  776275 
TOTAL  TEST  DAYS:  31267 
AVERAGE  MILES  PER  DAY:  24.83 
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APPENDIX  F 

WILCOX ON-RANK  SUM  STATISTICAL  ANALYSIS 
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¥ 


WI LCQX  ON-RANK  SUM  STATISTICAL  ANALYSIS 


The  ■first  test 

compared  populations  using  the  entire 

data  set.  Note:  The  two-tailed  probability  <P>  values  are 

compared  to  alpha  ■  .05/2  »  .025. 

TEST  ONE 

Factor:  Oil 

Type 

OIL  TYPE 

NUMBER 

2-TAILED  P 

STAUFFER 

37 

.2977 

MINERAL 

42 

STAUFFER 

37 

.1583 

CONOCO 

37 

MINERAL 

42 

.7162 

CONOCO 

37 

Factor:  Make  o-f  Vehicle 

MAKE 

NUMBER 

2-TAILED  P 

AMC 

14 

.1240 

CHEVROLET 

46 

AMC 

14 

.0622 

DODGE 

24 

AMC 

14 

.3989 

FORD 

17 

fitlC 

14 

.3426 

IHC 

5 

AMC 

14 

.9333 

OLDSMOBILE 

1 

f*1C 

14 

.7813 

PLYMOUTH 

9 

CHEVROLET 

46 

.3790 

DODGE 

24 

77 


CHEVROLET 

46 

.5986 

FORD 

17 

CHEVROLET 

46 

.6115 

IHC 

5 

CHEVROLET 

46 

.6383 

OLDSMOBILE 

1 

CHEVROLET 

46 

.1014 

PLYMOUTH 

9 

DODGE 

24 

.2040 

FORD 

17 

DODGE 

24 

.4476 

IHC 

5 

DODGE 

24 

.5600 

OLDSMOBILE 

1 

DODGE 

24 

.0785 

PLYMOUTH 

9 

FORD 

17 

.8795 

IHC 

5 

FORD 

17 

.8889 

OLDSMOBILE 

1 

FORD 

17 

.3666 

PLYMOUTH 

9 

IHC 

5 

.3333 

OLDSMOBILE 

1 

IHC 

5 

.2398 

PLYMOUTH 

9 

OLDSMOBILE 

1 

1 .0000 

PLYMOUTH 

9 

Factor:  Utilization  of  Veh 

i  cl  es 

USAGE  CODE 

NUMBER 

2-TAILED  P 

HI  UTE.LO  MI 

41 

.5279 

HI  UTE,HI  MI 

54 

HI  UTE,LO  MI 

41 

.0286 

LO  UTE,LO  MI 

18 
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HI 

UTEfLO 

MI 

41 

.1306 

LO 

UTE,HI 

MI 

3 

HI 

UTE.HI 

MI 

54 

.0097  * 

LO 

UTE.LO 

MI 

18 

HI 

LTTE, HI 

MI 

54 

.1560 

LO 

UTE.HI 

MI 

3 

LO 

UTE.LO 

MI 

18 

.0105  * 

LO 

UTE.HI 

MI 

3 

7? 


C  (low 


TEST  TUO 

Test  two  excludes  those  observations  with  Usage  code 
utilization/low  mileage). 


Factors  Oil 

type 

OIL  TYPE 

NUMBER 

2-TAILED 

STAUFFER 

30 

.4926 

MINERAL 

34 

STAUFFER 

30 

.2879 

CONOCO 

34 

, 

MINERAL 

34 

.7779 

CONOCO 

34 

Factor:  Make 

o-f  Vehicle 

MAKE 

NUMBER 

2-TAILED 

AMC 

11 

.0296 

CHEVROLET 

40 

AMC 

11 

.0590 

DODGE 

17 

AMC 

11 

.0712 

FORD 

16 

PMC 

11 

.0517 

IHC 

5 

AMC 

1 1 

.7103 

PLYMOUTH 

9 

CHEVROLET 

40 

.8753 

DODGE 

17 

CHEVROLET 

40 

.9855 

FORD 

16 

CHEVROLET 

40 

.9581 

IHC 

5 
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CHEVROLET 

40 

PLYMOUTH 

9 

0006E 

17 

FORD 

14 

DODGE 

17 

IHC 

3 

DODGE 

17 

PLYMOUTH 

9 

FORD 

14 

IHC 

5 

FORD 

14 

PLYMOUTH 

9 

IHC 

5 

PLYMOUTH 

9 

.  1-620 

.7898 

.9394 

.2409 

.7190 

.3014 

.2398 
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TEST  THREE 


Test  three  includes  only  those  observations  with 
Usage  code  C  (low  utilization/low  mileage). 


Factor]  Oil  Type 


OIL  TYPE 

NUMBER 

2-TAILED 

STAUFFER 

7 

.3357 

MINERAL 

8 

STAUFFER 

7 

1  .0000 

CONOCO 

3 

MINERAL 

8 

.4970 

CONOCO 

3 

Factor:  Make 

of  Vehicle 

MAKE 

NUMBER 

2-TAILED 

AMC 

3 

.7143 

CHEVROLET 

6 

AMC 

3 

.8333 

DODGE 

7 

AMC 

3 

.5000 

FORD 

1 

AMC 

3 

.5000 

OLDSMOBILE 

1 

CHEVROLET 

6 

.2343 

DODGE 

7 

CHEVROLET 

6 

.2857 

FORD 

1 

CHEVROLET 

6 

.2857 

OLDSMOBILE 

1 

DODGE 

7 

.2500 

FORD 

1 
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DODGE 

OLDSMOBILE 

FORD 

OLDSMOBILE 


7  .2500 

1 


1  1.0000 
1 
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APPENDIX  6 
fortran  program 


non  ooooonoonnonoooooooooooooooooooooooo 


PROGRAM  BLIND 


APPLICATION  OF  THE  TALBOTT  MODEL 
PEARCE  AND  DOUMIT;  SEP  1983 

SYMBOL  TABLE 


REAL: 

FAIL  -  OBSERVED  LIFETIMES  < UNORDERED) 

LIFE(I)  -  OBSERVED  LIFETIMES  (ORDERED) 

T< I )  -  CUMULATIVE  LIFE  AT  POINT  IN  TIME  LIFE(I) 

TAN  -  VALUE  OF  LIFE(I)  DEFINED  BY  POINT  OF  TANGENCY 
(THE  OPTIMAL  INTERVAL) 

U(I>  -  TTT  SCALED  STATISTIC;  T(I)/T(N);  Y-AXIS 
VARIABLE 

X(I)  -  LI FE( I )/LI FE(N) ;  X-AXIS  VARIABLE 
SL( I )  -  SLOPE  OF  THE  LINE  CONNECTING  K/LI FE(N)  TO 
THE  TTT  DATA  POINT 
MSL  -  MAXIMUM  VALUE  OF  SL(I) 

K  -  TIME  TO  REPLACE 
C  -  CRITICAL  TTT  POINT:  K/LIFE(N) 

B,E,Q,V,UJ,Y,Z  ,WW,ZZ  -  RANDOM  REAL  VALUES 

INTEGER 

I  -  OBSERVAT I ON  NUMBER 
N , J  —  INDICES 


FILES 

DATA  -  INPUT  DATA  FILE  CONTAINING  OBSERVED  OIL  LI  FES 
AND  OTHER  DESCRIPTIVE  DATA 
TRACK  -  OUTPUT  REPORT  FILE 


INITIALIZE  PROGRAM 

REAL  FAIL, LIFE( 1 50 ) ,T( 1 50 ) , TAN, U(150) , 
%  X(150> ,SL(150) ,MSL,K,C,W,WU,Y, 

%  ZZ , Z , B , E , 0 , V 
INTEGER  I,N,J 
OPEN  (10 , FILE—' DATA' ) 

OPEN  (11  ,  FI  LE58 'TRACK'  ) 

REWIND  10 
REWIND  11 


READ  IN  VALUES  FOR  N  AND  K 
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oooo  oooo  nnoo  oono 


READ< 1 0 , *)N , K 
WRITE  (11,800) 
WRITE  <11 , 820 )N , K 
WRITE  (11,810) 


READ  IN  <*4D  ORDER  LIFETIMES 

READd0,*)0,B,Z,W,Y,O,E,WU,ZZ,LIFE< 1) 
DO  20  1=2, N 

READ(10,*)0,B,Z ,W,Y,0,E,ZZ , YY ,FAI L 
LI FE< I )=FAI L 
DO  10  J= 1,1-1 

IF  (LIFE(I-J)  .ST.  FAIL)  THEN 
LIFE<1-J+1  )=LIFEd-J) 
LIFE(I-J)=FAIL 
ELSE 

GO  TO  20 
END  IF 
10  CONTINUE 

20  CONTINUE 

C=K/LIFE(N) 


COMPOTE  TTT  STATISTIC 

T< 1 )=LI FE< 1 ) *N 
DO  30  1*2, N 

T(I)=((LIFE<I) ~LI FE< I -1 > ) *<N-I + 1 .))+T(I-l) 
30  CONTINUE 


COMPOTE  AXIS  VARIABLES  AND  FIND  MAX  SLOPE 

MSL=Q . 

DO  40  1=1 ,N 

U< I )=T  < I )/T  <N) 

X( I )=LI FE< I )/LIFE(N) 

SL< I )=U< I )/<  C+X( I > ) 

WRITEdl  ,920)1  ,LIFE(I)  ,T(I)  ,U(I)  ,X(I)  ,SL<1) 
IF  < SL< I >  .GT.  MSL)  THEN 
MSL=SL(I) 

TAI^LI  FE<  I  > 

END  IF 
40  CONTINUE 

WRITEdl  ,?30)TAN 


FORMAT  STATEMENTS 
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AD  -  A  1 34  470 
UNCLASSIFIED 


OPTIMAL  MAINTENANCE  POLICIES:  A  GRAPHICAL  ANALYSIS(U)  4/'% 
AIR  FORCE  INST  OF  TECH  WRIGHT  PATTERSON  AFB  OH  SCHOOL  * 

OF  SYSTEMS  AND  LOGISTICS  P  F  DOUMIT  ET  AL .  SEP  83 
AFIT-LSSR-28-83  F/G  12/1 


NL 


800  FORMAT<///T20,  'TOTAL  TIME  ON  TEST'//) 

810  FORMAT <//T3,'I'  ,T8,'LIFE'  ,T19,  'TI '  ,T28,  'UI  '  , 

0  T34,'X<I>'  #T43,'SL'/> 

820  FORMATS  N-',I3,'  K -',F4.2> 

920  F0RMAT<T2I I3,T7f F4.0 ,T1?,F7.0 1T24,F3.3,T31 ,F7.3, 

♦  T40.F8.4) 

930  FORMAT <//' OPTIMAL  INTERVAL  IS  ',F4.0,'  MILES') 

STOP 
END 
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APPENDIX  H 

TTT  RESULTS i  ALL  OILS  COMBINED 
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TOTAL  TIME  ON  TEST 
ALL  OILS  COMBINED 


l*>116  K«  24.83 


I 

LIFE< I > 

T<  I  > 

U(I> 

X<I> 

SLOPE 

1 

707. 

82012. 

.104 

.041 

2.511747 

2 

1139. 

131492. 

.170 

.045 

2.534171 

3 

1971  . 

224540 . 

.292 

.113 

2.544074  * 

4 

2039. 

234224. 

.302 

.117 

2.543702 

3 

2599. 

294944. 

.383 

.149 

2.534574 

4 

2424. 

299719. 

.384 

.151 

2.534115 

7 

2894. 

329439. 

.425 

.144 

2.529537 

8 

3093. 

351112. 

.452 

.178 

2.524073 

9 

3284. 

371740. 

.479 

.189 

2.518103 

10 

3322. 

375804. 

.484 

.1 91 

2.514742 

11 

3330. 

374454. 

.485 

.191 

2.514404 

12 

3334. 

377074. 

.484 

.192 

2.514212 

13 

3481  . 

392342. 

.505 

.200 

2.508444 

14 

3587. 

403280. 

.520 

.204 

2.502581 

13 

3442. 

408890 . 

.527 

.209 

2.499335 

14 

3717. 

414445. 

.534 

.214 

2.494413 

17 

3720. 

414745. 

.537 

.214 

2.494410 

18 

3744. 

419141 . 

.540 

.215 

2.492454 

19 

3857. 

430215. 

.554 

.222 

2.484035 

20 

3914. 

435938. 

.542 

.225 

2.479395 

21 

3942. 

440354. 

.547 

.228 

2.475414 

22 

3972. 

441304. 

.548 

.228 

2.474747 

23 

3984. 

442432. 

.570 

.229 

2.473444 

24 

3988. 

442804. 

.570 

.229 

2.473258 

25 

3994. 

443540 . 

.571 

.230 

2.472440 

26 

3999. 

443813. 

.572 

.230 

2.472117 

27 

4022. 

445883. 

.574 

.231 

2.449532 

28 

4054. 

448731 . 

.578 

.233 

2.445807 

29 

4074. 

450447. 

.581 

.234 

2.443140 

30 

4138. 

454041 . 

.587 

.238 

2.455517 

31 

4303. 

470251 . 

.404 

.247 

2.435388 

32 

4442. 

499044. 

.443 

.247 

2.394871 

33 

4778. 

510490. 

.458 

.275 

2.382312 

34 

4874. 

518458. 

.448 

.280 

2.372083 

35 

4949. 

524408. 

.474 

.284 

2.344028 

34 

4970. 

524309. 

.478 

.284 

2.341722 

37 

4975. 

524709. 

.479 

.284 

2.341153 

38 

5083. 

535241 . 

.489 

.292 

2.348448 

39 

5104. 

534879. 

.492 

.293 

2.344210 

40 

5132. 

539035. 

.494 

.295 

2.342841 

41 

5154. 

540859. 

.497 

.294 

2.339879 

42 

5284. 

550459. 

.709 

.304 

2.323993 

43 

5298. 

551495. 

.710 

.305 

2.322243 
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44 

5323. 

553320 

45 

5326. 

553536 

44 

5341 . 

55460 1 

47 

5409. 

559361 

48 

5440. 

561500 

49 

5443. 

561704 

50 

5534. 

567801 

51 

5551 . 

548923 

52 

5602. 

572238 

53 

5632. 

574158 

54 

5728. 

580206 

55 

5782. 

583554 

56 

5887. 

589959 

57 

5947. 

593559 

58 

5967. 

594739 

59 

5991. 

596131 

60 

6060. 

600064 

41 

6085. 

601444 

62 

6104. 

602509 

63 

6215. 

608503 

64 

6225. 

609033 

65 

6273. 

611529 

66 

4379. 

616935 

47 

6414. 

418785 

68 

6457. 

420794. 

69 

6529. 

624250 

70 

6583. 

626788. 

71 

6603. 

627708. 

72 

6823. 

637608. 

73 

6964. 

643812. 

74 

6973. 

644199. 

75 

7021  . 

646215. 

74 

7153. 

651627. 

77 

7239. 

655067. 

78 

7351. 

659435. 

79 

7381. 

660575. 

80 

7435. 

662573. 

81 

7479. 

664157. 

82 

7579. 

667657. 

83 

7974. 

681087. 

84 

8105. 

685410. 

85 

8191  . 

688162. 

86 

8238. 

689619. 

87 

8312. 

691839. 

88 

8428. 

695203. 

89 

8444. 

695651 . 

90 

8508. 

697379. 

91 

8513. 

697509. 

92 

8578. 

699134. 

93 

8695. 

701942. 

94 

8809. 

704564. 

.713 

.306 

2.319036 

.713 

.306 

2.318640 

.714 

.307 

2.316607 

.721 

.311 

2.307251 

.723 

.313 

2.302936 

.724 

.313 

2.302508 

.731 

.318 

2.289399 

.733 

.319 

2 . 286929 

.737 

.322 

2.279405 

.740 

.324 

2.274924 

.747 

.329 

2.260525 

.752 

.332 

2.252426 

.760 

.338 

2.236704 

.765 

.342 

2.227743 

.764 

.343 

2.224721 

.768 

.344 

2.221032 

.773 

.348 

2.210334 

.775 

.350 

2.206425 

.776 

.351 

2.203407 

.784 

.357 

2.185741 

.785 

.358 

2.184144 

.788 

.361 

2.176380 

.795 

.367 

2.159277 

.797 

.369 

2.153310 

.800 

.371 

2.146637 

.804 

.375 

2.134873 

.807 

.378 

2.126035 

.809 

.380 

2.122731 

821 

.392 

2.086938 

.829 

.400 

2.064730 

830 

.401 

2.063314 

832 

.404 

2.055671 

839 

.411 

2.034767 

844 

.416 

2.021291 

849 

.422 

2.003871 

851 

.424 

1 .999204 

854 

.427 

1 .990735 

856 

.430 

1  .983793 

860 

.436 

1  .968021 

8  77 

.458 

1  .908468 

883 

.466 

1 .889634 

886 

.471 

1 .877362 

888 

.473 

1 .870635 

891 

.478 

1 .859999 

896 

.484 

1 .843394 

896 

.485 

1 .841097 

898 

.489 

1 .831827 

899 

.489 

1 .831096 

901 

.493 

1 .821494 

904 

.500 

1 .804272 

908 

.506 

1 .787640 
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95 

8844. 

705334. 

.909 

.508 

1 .782531 

94 

8854. 

705584. 

.909 

.509 

1.780759 

97 

8848. 

705824. 

.909 

.510 

1 .778941 

98 

8915. 

704719. 

.910 

.512 

1.771847 

99 

9795. 

722559. 

.931 

.543 

1.449218 

100 

10155. 

728479. 

.939 

.584 

1.404370 

101 

10248. 

730487. 

.941 

.590 

1.590493 

102 

10382. 

732197. 

.943 

.597 

1 .574951 

103 

10432. 

732897. 

.944 

.400 

1.570911 

104 

10821 . 

737954 . 

.951 

.422 

1 .525019 

105 

11245. 

743042. 

.957 

.444 

1 .477743 

104 

11494. 

747981 . 

.944 

.472 

1 .430589 

107 

11771. 

748751 . 

.945 

.477 

1 .422714 

108 

12384. 

754284. 

.972 

.712 

1 .342210 

109 

12800. 

757598. 

.974 

.734 

1 .324024 

110 

12931. 

758515. 

.977 

.743 

1 .312223 

111 

14884. 

770245. 

.992 

.854 

1.157804 

112 

15000. 

770815. 

.993 

.842 

1.149871 

113 

15512. 

772843. 

.994 

.892 

1.114933 

114 

15720. 

773487. 

.994 

.904 

1.101092 

115 

14829. 

775705. 

.999 

.947 

1 .031589 

114 

17399. 

774275. 

1  .000 

1 .000 

.998575 

OPTirttL  INTERVAL  IS  1971.  MILES 
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TOTAL  TIME  ON  TEST  PLOTS 
ALL  OILS  COMBINED 
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APPENDIX  I 

1 


TTT  RESULTS I  STAUFFER  SYNTHETIC  OIL 


TOTAL  TIME  ON  TEST 
STAUFFER  SYNTHETIC  OIL 


N-  37  K«  20.92 


I 

LIFE(  I ) 

T<  I ) 

U<I> 

X(I> 

SLOPE 

1 

707. 

24199. 

.119 

.042 

2.444454 

2 

1139. 

41711 . 

.184 

.048 

2.444203  * 

3 

1971. 

70831 . 

.312 

.117 

2.434444 

4 

2999. 

92183. 

.404 

.194 

2.408794 

S 

2894. 

101984. 

.449 

.172 

2.592258 

4 

3093. 

108288. 

.477 

.184 

2.978340 

7 

3284. 

114209. 

.903 

.195 

2.942183 

8 

3744. 

128009. 

.944 

.222 

2.520899 

9 

3942. 

134331. 

.992 

.235 

2.500407 

10 

3984. 

134947. 

.999 

.237 

2.498274 

11 

3994. 

139271. 

.994 

.237 

2.494791 

12 

4022. 

139947. 

.999 

.239 

2.493132 

13 

4074. 

137297. 

.405 

.242 

2.484702 

14 

4949. 

198249. 

.497 

.294 

2.340817 

19 

4970. 

198732. 

.499 

.295 

2.358059 

14 

9284. 

1 49440 . 

.730 

.314 

2.315033 

17 

9323. 

144499. 

.733 

.314 

2.309500 

18 

9934. 

170479. 

.752 

.329 

2.278101 

19 

9991  . 

171002. 

.753 

.330 

2.275449 

20 

9432. 

172440. 

.740 

.335 

2.241947 

21 

9947. 

178199. 

.785 

.355 

2.205935 

22 

4040. 

179443. 

.791 

.340 

2.190341 

23 

4219. 

181948. 

.802 

.349 

2.143541 

24 

4229. 

182108. 

.802 

.370 

2.141739 

29 

4379. 

184110. 

.811 

.379 

2.132915 

24 

4973. 

191238. 

.843 

.414 

2.027329 

27 

7021. 

191744. 

.845 

.417 

2.019049 

28 

7193. 

193084. 

.891 

.425 

1 .995541 

29 

7974. 

200479. 

.883 

.474 

1 .859140 

30 

8913. 

204787. 

.902 

.504 

1 .779199 

31 

8809. 

204899. 

.911 

.523 

1 .734954 

32 

8848. 

207213. 

.913 

.527 

1 . 728378 

33 

8919. 

207448. 

.914 

.530 

1 .721237 

34 

10432. 

213914. 

.941 

.420 

1 .514480 

39 

12931 . 

221013. 

.974 

.748 

1 .245184 

34 

19000. 

229191 . 

.992 

.891 

1.111 343 

37 

14829. 

224980 . 

1  .000 

1  .000 

.998758 

OPTIMAL  INTERVAL  IS  1139.  MILES 
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TOTAL  TIME  ON  TEST  PLOTS 
.  STAUFFER  SYNTHETIC  OIL 


APPENDIX  J 

TTT  RESULTS:  MINERAL  OIL 
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TOTAL  TIME  ON  TEST 
MINERAL  OIL 


N-  42  K>  28.30 


I 

LIFE(I) 

T<  I  > 

U(I> 

X<I> 

SLOPE 

1 

2039. 

85438. 

.292 

.130 

2.221549  * 

2 

3322. 

138241 . 

.472 

.21 1 

2.212839 

3 

3330. 

138541 . 

.473 

.212 

2.212478 

4 

3334. 

138717. 

.473 

.212 

2.212534 

S 

3442. 

150421 . 

.513 

.232 

2.197878 

4 

3717. 

153194. 

.523 

.234 

2.193400 

7 

3720. 

153304. 

.523 

.237 

2.193390 

8 

3857. 

158099. 

.539 

.245 

2 . 1 82234 

9 

4138. 

147453. 

.572 

.243 

2.158030 

10 

4442. 

184285. 

.429 

.295 

2.114128 

11 

4778. 

188437. 

.444 

.304 

2.104809 

12 

4874. 

191413. 

.454 

.310 

2.094147 

13 

5083. 

197883. 

.475 

.323 

2.074222 

14 

5154. 

200000. 

.482 

.328 

2.048884 

IS 

5298. 

203974. 

.494 

.337 

2.053742 

14 

5341  . 

205137. 

.700 

.340 

2.048911 

17 

5409. 

204905. 

.704 

.344 

2.040725 

18 

5440. 

207480 . 

.708 

.344 

2.034754 

19 

5443. 

207752. 

.709 

.344 

2.034345 

20 

5402. 

211409. 

.721 

.354 

2.013472 

21 

5782. 

215349. 

.735 

.348 

1 .987840 

22 

5887. 

217574. 

.742 

.374 

1 .972545 

23 

5947. 

218774. 

.744 

.378 

1 .943508 

24 

4414. 

227485. 

.777 

.408 

1 .894745 

23 

4403. 

231051 . 

.788 

.420 

1 . 848555 

24 

7239. 

241843. 

.825 

.440 

1 .784814 

27 

7351  . 

243455. 

.831 

.448 

1  .770748 

28 

7381  . 

244105. 

.833 

.470 

1 . 744834 

29 

7579. 

244877 . 

.842 

.482 

1  .740391 

30 

8105. 

253715. 

.844 

.514 

1  .472923 

31 

8191  . 

254747. 

.849 

.521 

1.442153 

32 

8428. 

257354 . 

.878 

.534 

1 .432102 

33 

8578. 

258854. 

.883 

.544 

1  .413003 

34 

8854. 

241354. 

.892 

.543 

1  . 577433 

35 

10248. 

272452. 

.930 

.453 

1 .420117 

34 

10382. 

273450 . 

.933 

.440 

1 .408477 

37 

11494. 

281322. 

.940 

.744 

1 .287027 

38 

11771  . 

281707. 

.941 

.749 

1 .280378 

39 

12384. 

284147. 

.949 

.788 

1  .227575 

40 

14884. 

291447. 

.995 

.947 

1 .048772 

97 


41  15512.  292919.  .999  .987 

42  15720.  293127.  1.000  1.000 

OPTIMAL  IWTERWAL  IS  2039.  MILES 


1 .010844 
.998203 
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TOTAL  TIME  ON  TEST  PLOTS 
MINERAL  OIL 


APPENDIX  K 

TTT  RESULTS*  CONOCO  SYNTHETIC  OIL 


100 


TOTAL  TIME  ON  TEST 
CONOCO  SYNTHETIC  OIL 


N-  37  K“  25.46 


1 

LI FE( I  ) 

T<  I  > 

U<I> 

X<I> 

SLOPE 

1 

2624. 

97088. 

.379 

.151 

2.488901  * 

2 

3481  . 

1 27940 . 

.499 

.200 

2.478203 

3 

3587. 

131650. 

.514 

.206 

2.475240 

4 

3916. 

142836. 

.558 

.225 

2.461387 

5 

3972. 

144684. 

.565 

.228 

2.458305 

6 

3988. 

145196. 

.567 

.229 

2.457170 

7 

3999. 

145537. 

.568 

.230 

2.456208 

8 

4054. 

147187. 

.575 

.233 

2.450565 

9 

4303. 

154408. 

.603 

.247 

2.422902 

10 

4975. 

173224. 

.676 

.286 

2 . 352868 

11 

5104. 

176707. 

.690 

.293 

2.339815 

12 

5132. 

177435. 

.693 

.295 

2.336699 

13 

5326. 

182285. 

.712 

.306 

2.313546 

14 

5728. 

191933. 

.749 

.329 

2.265792 

IS 

5991  . 

197982. 

.773 

.344 

2.235034 

16 

6085. 

200050. 

.781 

.350 

2.223638 

17 

6104. 

200449. 

.782 

.351 

2.221166 

18 

6273. 

203829. 

.796 

.361 

2.198017 

19 

6457. 

207325. 

.809 

.371 

2.172257 

20 

6529. 

208621 . 

.814 

.375 

2.161825 

21 

6583. 

209539. 

.818 

.378 

2.153595 

22 

6823. 

213379. 

.833 

.392 

2.116207 

23 

6964. 

215494. 

.841 

.400 

2.094069 

24 

7435. 

222088. 

.867 

.427 

2.021896 

25 

7479. 

222660 . 

.869 

.430 

2.015219 

26 

8238. 

231768. 

.905 

.473 

1 .904982 

27 

8312. 

232582. 

.908 

.478 

1 .894706 

28 

8444. 

233902. 

.913 

.485 

1 . 875762 

29 

8508. 

234478 . 

.915 

.489 

1 . 866278 

30 

8695. 

235974. 

.921 

.500 

1 .837910 

31 

8844. 

237017. 

.925 

.508 

1 .815021 

32 

9795. 

242723. 

.948 

.563 

1 .678721 

33 

10155. 

244523. 

.955 

.584 

1 .631367 

34 

10821 . 

247187. 

.965 

.622 

1 . 547879 

35 

11245. 

248459. 

.970 

.646 

1 .497313 

36 

12800. 

251569. 

.982 

.736 

1 . 332243 

37 

17399. 

256168. 

1.000 

1  .000 

.  998539 

OPTIMAL  INTERVAL  IS  2624.  MILES 
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